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ABSTRACT 

/  ^  I 

In  this  paper  we  study  the  steady-state  behavior  of  solids  that  can 
sustain  mechanical/  electromagnetic/  and  thermal  effects*  examine  a  class 
of  boundary-value  problems  for  a  guasilinear  system  of  functional  differential 
equations  that  is  derived  from  a  very  general  model  for  such  materials. 
propose  a  physically  reasonable  constitutive  theory  which  leaves  this  system 
amenable  to  modem  methods  of  partial  differential  equations.  The  principal 
assumption  is  a  modified  version  of  the  strong  elliptlcity  condition. ^liv^ Part 
I  we  prove^existence  results  for  the  general  system  under  some  special 
physical  assumptions  (rigidity  and  hyperelasticity) .  to  formulation  admits 
non-local  interactions  caused  by  the  magnetic^" self-field"igenerated  by  the 
deformed,  conducting  body.  ~£n^¥art  It  we  •fehow^the  existence  and  regularity  of 
solutions  of  a  system  of  functional  ordinary  differential  equations  arising 
from  a  semi-inverse  problem  in  a  more  comprehensive  physical  situation. 
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SIGNIFICANCE  AND  EXPLANATION 


In  this  paper  we  study  the  steady-state  behavior  of  solids  that  can 
sustain  mechanical,  electromagnetic,  and  thermal  effects.  Our  goal  is  to 
formulate  very  general  assumptions  on  the  constitutive  equations  that  are 
physically  reasonable,  yet  leave  the  resulting  mathematical  problems 
tractable.  The  models  we  propose  admit  the  following  types  of  nonlinear 
behavior  which  are  of  particular  interest  when  materials  are  subjected  to 
large  electromagnetic  fields  and  sustain  large  currents. 

1.  Nonlinear  Coupling:  We  investigate  the  abstract  mathematical  problems 

that  occur  when  we  assume  a  very  general  coupling  of  the  various  physical 
fields  of  the  material.  For  instance,  we  assume  that  the  dielectric 
displacement  d  depends  not  only  on  the  electric  field  $  but  also  on  the 
deformation  gradient  £,  temperature,  temperature  gradient  jg,  and  magnetic 
field  h.  We  make  similar  assumptions  about  the  stress  tensor  £,  magnetic 
induction  h,  heat  flux  vector  g,  and  electric  current  The  assumption 

which  makes  such  problems  tractable  is  a  modified  version  of  the  "strong 
ellipticity  condition". 

2.  Nonlocal  Self-Interactions i  The  Maxwell's  equation 

curl  h  -  1 

fSf 

Implies  that  a  current  in  one  part  of  a  body  will  generate  *  magnetic  field  in 
a  distant  part  of  the  body.  Since  we  allow  J  to  depend  on  the  entire  list 
of  independent  variables,  the  magnetic  field  at  any  point  will  depend  on  the 
global  values  of  the  other  variables.  And,  since  all  of  our  dependent 
variables  depend  on  h,  this  problem  will  be  spread  to  the  entire  system  of 
equations.  We  use  compactness  method  to  handle  these  nonlocal  problems. 

In  Part  I  we  handle  very  general  mathematical  problems  under  special 
physical  assumptions  (rigidity  and  hyperelasticity)  and  in  Part  II  we  handle  a 
special  mathematical  problem  (a  semi-inverse  problem)  under  general  physical 
conditions. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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STEADY -STATE  PROBLEMS  OF  NONLINEAR  ELECTRO-MAGNETO-THERMO- ELASTICITY* 

1  2 
Robert  C.  Rogers  and  Stuart  S.  Antman 

Part  I.  The  General  Theory 

1 .  Introduction 

In  this  paper  we  study  a  class  of  boundary  value  problems  for  a  quasi linear 
system  of  functional-differential  equations  describing  the  steady-state  behavior 
of  solids  that  can  sustain  mechanical,  electromagnetic,  and  thermal  effects.  We 
treat  partial  differential  equations  in  Part  I  and  ordinary  differential 
equations  in  Part  II.  Our  primary  goals  are  to  show  that  there  is  a  simple  way 
to  formulate  the  governing  equations,  which  illuminates  the  physics  and  promotes 
the  analysis  of  the  equations,  to  actually  analyze  important  classes  of  problems, 
and  to  contribute  to  the  development  of  an  effective  constitutive  theory  for  such 
materials  by  showing  how  our  physically  and  mathematically  natural  constitutive 
restrictions  support  existence  and  regularity  theories  for  our  problems.  The 
problems  we  study  are  simple  enough  to  be  tractable,  interesting  enough  to 
possess  a  very  rich  class  of  solutions,  and  yet  complicated  enough  to  require  new 
approaches,  both  in  the  formulation  and  treatment  of  electromagnetism  in  solids 
and  in  the  use  of  techniques  of  nonlinear  analysis. 

Our  constitutive  equations  give  the  stress,  heat  flux,  dielectric 
displacement,  magnetic  induction,  and  electric  current  as  arbitrary  functions  of 

*Some  of  the  results  reported  here  were  developed  in  the  doctoral  dissertation  of 
Rogers  (1984). 

^Partially  sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C- 
0041.  This  material  is  partially  based  upon  work  supported  by  the  National 
Science  Foundation  under  Grant  No.  DMS-82 10950,  Mod.  1. 

^partially  supported  by  National  Science  Foundation  Grant  No.  DMS-8503317. 


the  deformation  gradient,  temperature,  temperature  gradient,  electric  field,  and 
magnetic  field.  These  constitutive  functions  must  of  course  be  invariant  under 
rigid  motions.  In  order  to  reduce  the  governing  equation  to  ordinary 
differential  equations  for  our  semi-inverse  problems  of  Part  II,  we  further 
require  that  the  bodies  under  study  have  some  material  symmetry.  Our  basic 
constitutive  assumption  is  that  the  constitutive  equations  satisfy  the  Strong 
Ellipticity  Condition.  Originating  in  the  theory  of  partial  differential 
equations,  this  condition  proves  to  be  eminently  natural  on  physical  grounds. 
Indeed,  this  condition,  when  precisely  formulated,  is  equivalent  to  the 
requirement  that  each  component  of  the  dependent  constitutive  variables  is  a 
strictly  increasing  function  of  the  corresponding  component  of  the  independent 
constitutive  variables  (when  the  other  components  of  the  independent  variables 
are  held  fixed) .  Roughly  speaking,  a  typical  consequence  of  this  assumption  is 
that  a  change  in  the  temperature  gradient  produces  a  far  more  pronounced  change 
in  the  heat  flux  vector  than  it  does  in  the  stress,  dielectric  displacement,  and 
magnetic  induction.  Thus  the  Strong  Ellipticity  Condition  implies  a  very  mild 
uncoupling  in  the  constitutive  equations.  True  uncouplings,  such  as  the 
independence  of  stress,  dielectric  displacement,  and  magnetic  induction  on  the 
temperature  gradient,  may  be  interpreted  as  consequences  of  the  Clausius-Duhem 
version  of  the  Second  Law  of  Thermodynamics.  But  we  have  no  need  for  such  true 
uncouplings  anywhere  in  our  analysis.  Indeed,  wherever  the  Clausius-Duhem 
inequality  is  more  restrictive  than  the  Strong  Ellipticity  Condition,  we  have  no 
need  for  its  consequences,  and  wherever  it  is  less  restrictive,  it  is  inadequate 
for  our  needs.  We  supplement  the  Strong  Ellipticity  Condition  with  compatible 
growth  conditions. 

Rather  than  adhering  to  the  classical  tradition,  (followed  at  least  in  part 
by  Toupin  (1956),  Fano,  Chu,  &  Adler  (I960),  Penfield  &  Haus  (1967),  and  DeGroot 


&  Suttorp  (1972)  and  others ),  of  deriving  or  motivating  the  constitutive 
equations  of  electromagnetism  from  discrete  microscopic  models ,  we  employ  the 
phenomenological  approach  of  continuum  mechanics  and  simply  lay  down  general 
constitutive  laws.  We  thereby  gain  great  economy  and  generality  in  our 
formulation  of  electromagnetism.  (Cf.  Truesdell  &  Toupin  (I960).) 

For  simplicity  and  clarity  in  our  mathematical  analysis ,  it  is  crucial  not 
only  that  we  give  Maxwell's  equations  a  material  (Lagrangian)  formulation,  but 
also  that  we  introduce  new  variables  in  place  of  the  dielectric  displacement, 
magnetic  induction,  etc.  Zn  this  regard  we  are  merely  extending  to  the  theory  of 
electromagnetism  in  deformable  media  the  methodology  that  has  proved  most  natural 
and  successful  for  boundary  value  problems  of  nonlinear  elasticity. 

In  the  overwhelming  majority  of  texts  on  continuum  mechanics  the  emphasis 
placed  on  the  spatial  formulation  of  the  governing  equations  and,  in  particular, 
on  the  use  of  the  Cauchy  stress  tensor  overshadows  that  placed  on  the  material 
formulation  and  on  the  use  of  the  Piola-Kirchhoff  stress  tensors.  (We  define 
these  stress  tensors  in  Section  3.)  The  reason  for  this  emphasis  is  largely 
historical:  The  two  most  highly  developed  branches  of  continuum  mechanics  are 
Newtonian  fluids  and  linear  elasticity.  The  constitutive  equations  for  a 
homogeneous  Newtonian  fluid  are  especially  simple  in  a  spatial  formulation. 
Moreover,  in  this  formulation  the  constraint  of  incompressibility  (valid  for 
liquids)  has  an  elegant  characterization  as  the  linear  equation  expressing  the 
vanishing  of  the  divergence  of  the  velocity  field  defined  over  points  in  space. 
(In  contrast,  incompressibility  is  characterized  in  a  material  formulation  by  the 
nonlinear  equation  requiring  the  Jacobian  of  the  deformation  gradient  to  equal 
unity.)  For  problems  involving  nonhomogeneous  fluids  or  fluids  with  free 
surfaces,  there  are  compensating  disadvantages  requiring  some  version  of  a 
material  formulation,  possibly  disguised,  for  their  successful  treatment.  In 
linear  elasticity  there  is  no  distinction  between  material  and  spatial 
formulations,  although  the  derivation  of  these  equations  from  a  nonlinear  spatial 
formulation  is  much  more  difficult  to  carry  out  than  that  from  a  nonlinear 
material  formulation.  The  advantages  of  a  material  formulation  are  evident  for 
the  boundary  value  problems  of  nonlinear  solid  mechanics:  (i)  The  prescription 
of  constitutive  equation  for  the  first  Piola-Kirchhoff  stress  tensor  in  terms  of 
the  past  history  of  deformation  is  natural  and  does  not  suffer  from  complications 
due  to  nonhomogeneity,  (ii)  The  governing  equations  are  posed  on  a  fixed  region 
of  space,  the  region  occupied  by  the  body  in  a  reference  configuration,  rather 
than  on  the  unknown  and  possibly  moving  region  actually  occupied  by  the  body. 
These  factors  have  not,  however,  proved  to  be  compelling  in  shifting  the  emphasis 
of  texts  towards  material  formulations  because  there  have  been  so  few  studies  of 
nonlinear  boundary  value  problems  of  solid  mechanics.  (Cf.  Antman  (1978,  1979, 
1983),  Ball  (1977,  1982).) 


Maxwell's  equations  have  been  posed  almost  exclusively  in  spatial 
coordinates  because  the  most  important  case  of  the  vacuum  can  be  posed  in  no 
other  way  and  because  in  the  second  most  important  case  of  a  rigid  medium  there 
is  no  essential  distinction  between  material  and  spatial  coordinates.  Moreover, 
the  most  actively  cultivated  field  of  electromagnetism  in  deformable  media  is 
that  of  magnetohydrodynamics .  For  the  reasons  mentioned  in  our  connents  on  fluid 
dynamics,  many  problems  for  this  theory  are  most  easily  set  in  spatial 
coordinates.  Only  recently  has  the  use  of  material  coordinates  begun  to  appear 
in  treatments  of  electromagnetism  in  media.  (Cf.  Walker,  Pipkin,  6  Rivlin 
(1965),  Hutter  (1975),  Pao  &  Hutter  (1975),  Hutter  &  van  de  Ven  (1978),  McCarthy 
&  Tiersten  (1978),  and  Maugin  (1981).)  These  authors  have  also  introduced  new 
fields  suitable  for  material  coordinates  in  place  of  the  classical  fields. 

The  large  deformation  of  solids  in  the  presence  of  large  electromagnetic 
fields  is  a  problem  of  growing  technological  importance  (cf.  Moon  (1978,  1984)). 
Awareness  of  this  importance  is  evidenced  by  the  number  of  papers  recently 
devoted  to  this  subject.  (Cf.  Parkus  (1979),  Ambartsumian  (1982),  Maugin 
(1983).)  To  our  knowledge,  ours  is  the  first  mathematical  analysis  of  general 
nonlinear  boundary  value  problems  in  this  area. 

Much  of  the  previous  work  in  the  electromagnetism  of  deformable  solids  can 
be  divided  into  two  general  areas:  the  development  of  a  general  theory  governing 
such  media  and  the  solution  of  specific  nonlinear  problems.  General  theories  of 
the  dynamics  of  deformable  solids  have  been  proposed  by  Fano,  Chu  &  Adler  (1960), 
Toupin  (1963),  Dixon  &  Eringen  (1965),  Pao  &  Hutter  (1975),  and  Maugin  &  Eringen 
(1977).  These  developments  consist  in  the  derivation  of  some  form  of  Maxwell's 
equation  and  associated  forms  of  the  electric  body  force,  body  couple,  and 
internal  energy  supply  from  some  discrete  model  of  the  material.  Comparisons  of 
various  theories  are  to  be  found  in  Penfield  &  Haus  (1967),  DeGroot  &  Suttorp 
(1972),  Hutter  &  van  de  Ven  (1978),  and  Pao  (1978).  Our  work  is  more  closely 
related  to  the  static  theories  of  Toupin  (1956)  and  Brown  (1966).  Specific 
problems  for  general  nonlinear  dielectrics  were  solved  by  Toupin  (1956),  Eringen 
(1963),  Verma  (1964),  Pipkin  &  Rivlin  (1960),  and  Singh  &  Pipkin  (1966)  via  the 
inverse  methods  of  modem  nonlinear  elasticity.  (Singh  &  Pipkin  also  provide  a 
review  of  the  earlier  work.)  There  is  an  extensive  literature  on  specific 
nonlinear  materials  with  polynomial  constitutive  equations  and  associated 
problems  (cf.  Jordan  &  Eringen  (1964)  and  Pipkin  &  Rivlin  (1966)).  Maugin  (1981) 
reviews  the  modern  work  on  wave  motion  in  magnetizable  deformable  solids  and 
includes  both  general  nonlinear  and  specific  (approximate)  constitutive 
equations. 


2 .  Notation 


The  euclidean  3-apace  X3  is  defined  to  be  abstract  3-dimensional  real 
inner-product  space.  Its  inner  product,  the  dot  product,  is  the  natural  source 
of  the  geometric  properties  of  the  space.  We  interpret  X3  as  physical  space. 

We  distinguish  X3  from  X3,  the  space  of  triples  of  real  numbers  equipped  with 
any  norm  (which  is  necessarily  topologically  equivalent  to  the  Euclidean  norm) . 
But  we  assign  no  natural  geometrical  significance  to  the  norm  on  X  • 

Vectors ,  which  we  define  to  be  elements  of  X3,  and  vector-valued  functions 
are  denoted  by  bold-face,  lower-case  Latin  letters.  Second-order  tensors,  which 
form  the  space  Lin  of  linear  operators  from  X3  into  itself,  are  denoted  by 
bold-face  upper-case  Latin  letters.  The  subspace  of  Lin  consisting  of 
symmetric  second-order  tensors  is  denoted  Sym.  Its  subset  of  positive-definite 
tensors  is  denoted  Psym.  The  group  SL(3)  of  all  members  of  Lin  with 
positive  determinant  is  denoted  Lin+.  Scalars  and  scalar-valued  functions  are 
denoted  with  light-faced  letters.  Elements  of  for  n  >  1  and  functions 

with  values  in  r"  are  denoted  by  bold-face  sanserif  lower-case  Latin  letters 
and  by  bold-face  lower-case  Greek  letters. 

We  employ  the  dyadic  notation  of  Gibbs  (cf.  Gibbs  &  Wilson  (1901)),  which  we 
now  describe.  (This  notation  is  both  admirably  suited  for  treatment  of  problems 
in  curvilinear  coordinates  and  completely  compatible  with  modern  invariant 
formulations  of  linear  algebra  in  X3. )  The  dot  product  of  vectors  u  and  v 
is  denoted  The  cross  product  of  two  vectors  and  y  is  denoted  jjvy. 


The  value  of  the  second  order  tensor  at  y  is  denoted  &*sj*  The  transpose 

#  * 

A  of  ft  is  defined  by  v*(A*u)  *  u*(A  »v)  for  all  jj  and  We 

•  * 

accordingly  write  A  »v  *  v*A.  ^  is  symmetric  if  A  ■  A  and  skew  if 

A  »  -A  .  If  A  is  skew,  there  is  a  unique  vector  a,  called  the  axial  vector 


such  that 


for  all 


The  dyadic  product 


vectors  jj  and  v  is  defined  to  be  the  second-order  tensor  satisfying 

* 

<HX>*K  “  for  a11  58*  Thus  (uv)  -  vu  and  tr(jQj)  *  where  tr 

denotes  trace.  The  product  (ft»0)  of  tensors  is  defined  by  (ft*0)*X  ”  0*  (0*2) 

for  all  x.  Thus  -  (ft*jj)x  and  (»X)*(ft)  ■  }j(x*ft).  We  8et 

*  *  * 

A:B  ”  tr(A*B  ).  Hence  trA  ■  I:A  ■  A:I  ■  A  il  ■  trA  , 

A/  /V  *W  (V  f\J  /V  (V  fV 

(»X)S(XX)  *  (0*X) (X*X)*  and  A* (JUX)  “  W*A*X  ■  ()4X)8&.  (i  denotes  the  identity 
tensor.)  It  is  easy  to  see  that  is  an  inner  product  on  Lin.  We 

accordingly  define  the  norm  | A |  of  A  by  |a|  *  /AsA.  If  a  and  b  are  unit 
vectors,  then  |ab|  ■  1  so  that  ab  is  a  unit  tensor.  In  this  case  we  can 
represent  any  tensor  A  by  the  orthogonal  decomposition 
A  *  (a‘A‘b)ab  +  [A  -  ( a*A*b)gb]  where  a*  A*  b  is  the  component  of  A  along 
ab  and  [A  -  (a.ft*b)ab]  is  the  projection  of  A  onto  the  orthogonal 
complement  of  ab.  If  {jj^}  and  {jb^}  are  each  bases  for  K3,  then 
is  a  basis  for  Lin.  Thus  we  can  use  all  of  our  dyadic  identities  to  construct 
the  familiar  componential  formulas  for  all  the  expressions  we  have  introduced  in 
the  invariant  form.  Repeated  indices  are  summed  over  their  obvious  ranges. 

The  (Gateaux)  differentials  of  u  ^  f(u)  at  a  in  the  direction  of  £  and 
of  H  ►  0(5J)  at  ft  in  the  direction  of  0  are  defined  to  be  the  vector 
t3j£Ca>/3jj]  and  the  tensor  [30(ft)/3iJ]  *0  given  by 
(2.D  [3f(a)/3u).jb  -|^f(a  +  tb)|t_Q  , 

(2-2)  [3F(A)/30]«B  -|^F(A  +  tB)|t>t0  . 

Other  differentials  are  defined  similarly.  If  U  *►  F(U)  and  V  *►  G(V) 
are  (Frechet)  differentiable,  then  V  »  F(G(V))  «  H(V)  is  also,  and  its 
differential  satisfies  the  chain  rule: 

(2.3)  13H(A)/3V]:B  -  [3f(G( A) )/3u] : { [3G( A)/3V] :B>  . 

(The  braces  can  be  omitted  from  the  right  side  of  (2.3).)  Our  notational  scheme 


embodied  in  (2.1)  and  (2.2)  causes  the  chain  rule  (2.3)  to  have  a  form  analogous 


to  that  tor  scalar  functions*  As  we  shall  see  in  the  next  chapter  this  virtue  is 


counterbalanced  by  the  increased  complexity  of  defining  and  representing  the 
action  of  the  classical  differential  operators  grad,  div,  curl  on  tensor 
functions . 

We  denote  n  copies  of  a  function  space  X  by  X  itself.  The  distinction 
will  be  clear  from  the  contexts  Thus  a  statement  of  the  form  w  e  I.  (8)  is  to 

~  p 

imply  that  this  Lp(8)  is  the  space  of  all  measurable  vector-valued  functions 

D  8  3  z  *  w(z)  e  B3  such  that  /  [w(s)»w(z)]p,/,2dv(z)  <  •  . 

B 


The  norm  of  a  Banach  space  X  is  denoted  I*, XI 


3.  Formulation  of  the  Governing  Equations 


In  this  and  the  next  section  we  formulate  the  equations  for  steady-state 
problems  of  electro-magneto-thermo-elasticity.  There  are  several  different 
theories  that  are  at  least  formally  equivalent  in  the  classical  nonrelativistic 
setting  we  employ.  (Cf.  Hutter  &  van  de  Ven  (1978)  and  Pao  (1978).)  Comparisons 
of  the  various  theories  is  made  difficult  by  the  fact  that  the  same  symbol  used 
in  different  theories  has  different  meanings.  Fortunately,  the  mathematical  form 
of  the  governing  equations  expressing  the  balance  of  linear  momentum,  the  balance 
of  energy.  Maxwell's  laws,  and  the  conservation  of  charge  is  the  same  for  all 
these  theories.  We  shall  refer  to  the  various  fields  that  occur  in  our  equations 
by  their  traditional  names,  realizing  that  their  precise  physical  significance 
inheres  in  the  slots  they  occupy  in  the  equations  for  a  specific  theory. 

To  make  our  presentation  as  transparent  as  possible,  we  assume  that  all  the 
functions  and  boundaries  that  appear  are  smooth  enough  for  all  the  classical 
operations  that  appear  to  make  sense.  (A  careful  treatment  of  these  issues 
without  such  blanket  smoothness  assumptions  can  be  modelled  on  that  of  Antman  & 
Osborn  (1979).)  Of  course,  we  abandon  this  optimistic  formalism  when  we 
afterwards  analyze  our  specific  boundary  value  problems. 

3 

We  identify  a  material  body  with  the  closure  B  of  a  domain  in  E  and  we 
identify  material  points  of  the  body  with  their  positions  z  in  8.  For  each 
z  in  8  let  y(j)  denote  the  position  of  z  in  some  deformed 
configuration.  The  (transposed)  deformation  gradient  F  and  the  right  Cauchy- 
Green  deformation  tensor  £  for  the  configuration  y  are  defined  by 
(3.1)  F  =  3y/3z,  C  =  F* *F  . 

We  require  that  no  two  distinct  material  points  simultaneously  occupy  the 
same  position  in  a  given  configuration.  Thus  each  map  y  must  be  one-to-one. 
Since  this  global  condition  is  so  difficult  to  treat,  we  ignore  it  and  content 


ourselves  with  the  local  condition  that  the  deformation  £  merely  preserve 
orientation,  i.e.,  that 

(3.2)  det  F  >  0  , 
where  det  denotes  the  determinant. 

Let  X(y)  denote  the  logarithm  of  the  absolute  temperature  at  position  y 
in  space.  (It  is  finite-valued  if  and  only  if  the  absolute  temperature  is 
posit ive- valued. )  We  set  g(y)  *  3X(y)/3y.  Let  e(y)  and  h(y)  denote  the 
electric  and  magnetic  fields  at  y.  We  set 

(3.3)  X(S)  =  X(y(z)),  g(z)  =  3X(z)/3z  =  g(X(z))*F(z)  , 

S(£>  5>»  “  E<X<£>>*£<*>  • 

g,  e,  and  h  are  the  material  logarithmic  temperature  gradient,  electric 

field,  and  magnetic  field.  (We  shall  soon  see  that  e  and  h  can  be 
represented  in  terms  of  gradients.  Consequently  they  transform  in  (3.3)  just 
like  g. ) 

Let  T(y)  denote  the  effective  Cauchy  stress,  i.e.,  the  sum  of  the 
mechanical  Cauchy  stress  and  the  Maxwell  stress.  T(y)  measures  force  per  unit 
actual  area  at  y.  (There  are  several  versions  of  T,  depending  on  alternative 

A# 

representations  and  decompositions  of  the  Lorentz  force  and  torque.)  Let  g(y) 

a#  <w  *w 

denote  the  heat  flux  per  unit  actual  area  at  y.  Let  d(y),  b(y),  J(y)  be  the 
dielectric  displacement,  magnetic  induction,  and  current  density  at  y.  Then  we 
introduce  material  versions  of  these  fields  by 

(3.4)  T(z)*  -  det  F(z)F“1(z)»T(y(z))*  , 

g(z)  =*  det  F(z)F_1(z)*|(y(z)),  etc. 

T  is  the  effective  first  Piola-Kirchhoff  stress. 

- -  ■  -  -  -  -  • 

For  simplicity  let  us  assume  that  the  body  force  and  heat  source  have  purely 

a 

electromagnetic  origin.  Then  the  local  form  of  the  balance  of  forces,  the 


balance  of  energy,  and  Maxwell's  equations  for  a  steady  state  are 


(3.5) 


(3.6) 


(3.7) 

(3.8) 


(3.9) 


(3.10) 


Div  T  +  £  ■  0  , 

A»  (V  tv 

Div  3  +  4*S  "  0  < 
Div  £  -  a  , 
Div  b  -  0  , 

<v 

Curl  S  -  0  , 
Curl  fe  -  i  . 


The  material  divergence  Div  of  a  tensor  is  defined  by  Green's  Theorem 


(3.11) 


/  T*nda  ■  /  Div  T  dv 

'  tv  tv  *  A* 

3P  P 


where  PC  8  and  n  is  the  unit  outer  normal  to  P.  In  (3.5)  f  represents 
body  forces  of  electromagnetic  origin  not  absorbed  by  the  Maxwell  stress.  Since 
every  term  in  the  usual  prescriptions  of  the  Lorentz  force  is  a  divergence,  we 
could  absorb  this  force  entirely  into  the  Maxwell  stress  and  hence  into  the 
effective  stress.  We  accordingly  take  f  ■  0.  (Cf.  Hutter  6  van  de  Ven  (1978).) 


The  term  ^*e  in  (3.6)  is  the 


i.  In  (3.7)  or  represents  the  free 


charge.  We  regard  it  as  an  assigned  function  of  z. 


The  balance  of  torque  has  the  local  form 


(3.12) 


*  * 
L  -  T*F  -  F»T 


where  Jj  is  a  skew  tensor  depending  upon  the  electromagnetic  fields  and  the 
choice  of  the  Maxwell  stress  tensor.  We  assume  that  (3.12)  is  identically 


satisfied  when  the  constitutive  functions,  to  be  introduced  in  the  next  section, 
are  substituted  into  (3.12).  Hutter  &  van  de  Ven  (1978)  show  that  it  is 
permissible  to  take  L  ■  O. 

Equations  (3.9)  and  (3.10)  imply  that  there  exist  scalar  functions  9  and 
i|»,  called  the  electric  and  magnetic  scalar  potentials  such  that 


(3.13) 


e  -  89/8*  , 


V-V*V  V-V'-V-V*'.'  V 
4  4  .  4  . 


(3.14)  h(£)  -  3*(s>/3*  +  f  U(u)  -  (a  -  *)]|«  “  *r3<iv{2>  ' 


B 


as  xs 


shown  in  standard  books  on  electromagnetism.  (These  formulas  justify  the 


remarks  following  (3.3).) 


Constitutive 


Equations 

Of  all  the  variables  that  have  appeared  only  a  is  prescribed.  The 
remaining  variables  are  related  by  constitutive  equations.  As  our  independent 
constitutive  variables  we  choose 

(4.1)  £,  = 

because  they  are  physically  reasonable  and  mathematically  convenient.  The  domain 
of  F  is  Lin+,  the  domain  of  is  E3,  and  the  domain  of  X  is  R.  We 

first  suppose  that  <j,d,b,^  depend  on  these  variables  and  on  g.  Thus 
appearing  in  (3.6),  likewise  depends  on  (4.1).  We  finally  prescribe  T  and  £ 
to  depend  on  (4.1)  and  j.  so  that  they  satisfy  (3.12).  Henceforth  we  shall  have 
no  need  for  (3.12).  Thus  we  have  constitutive  functions  such  that 

(4.2)  T(z)  -  T(Hz),z),  etc. 

A 

The  functions  T,  etc.,  must  be  invariant  under  rigid  motions,  i.e.,  be  frame- 
indifferent.  (Cf.  Truesdell  &  Noll  (1965).)  We  do  not  pause  to  exhibit  the 
specific  representations  of  the  constitutive  functions  that  are  necessary  and 
sufficient  for  frame-indifference  because  we  have  no  need  for  them  in  our 
analysis. 

For  simplicity,  we  assume  that  our  constitutive  functions  are  continuously 


di f  f  erentiable 


ru  vn.  r.wura*'.  FlTiKK* 


5.  Potentials 

Xt  is  convenient  in  our  analysis  to  employ  the  potentials  p  and  $ 
instead  of  e  and  Ji  as  the  fundamental  variables  defining  the  electromagnetic 
state,  e  is  expressed  as  the  gradient  of  9  in  (5.15).  If  the  current 
^  -  0,  then  )rj  is  likewise  expressed  as  the  gradient  of  $  by  (5.16).  We 
seek  conditions  ensuring  that  jj  can  be  expressed  in  terms  of 


*s  »  *V 


(5.1) 


(JE'3A,3p/3;8,3^/3;g)  =  A 


when  the  current  is  not  sero.  Mote  that  each  entry  in  _A  except  X  is  a 
gradient.  Let  us  substitute  our  constitutive  equation  for  i  into  (3.14)  to 


obtain 

(5.2)  h(z)  -  3«(j£)/3£  “  /  ti(h(x),£(x),x)  -  (x  -  z)]|x  -  zr"3dv(x)  5  J$(h,£)(z> 

8 


w.v; 

\v.>: 

r.v.v] 
r  v  V 


where  l  =  (j;,jg,X,g).  Mow  in  the  classical  form  of  Ohm's  Law,  £  depends  only 

A 

on  the  electric  field.  More  generally,  if  ^  is  independent  of  h,  then  (5.2) 
gives  an  explicit  representation  for  in  terms  of  A*  *h«*e  are  a  variety  of 

A 

results  available  for  the  case  that  j  depends  on  h.  Typical  is  the  following: 
5.3.  Theorem.  Let  a  >  1  and  let  8  lie  in  the  ball  of  radius  Y  and 

center  g.  Let  A  be  fixed  in  La(8).  Suppose  that  there  are  positive  numbers 
y,  0,  C  with  3C  <  o  such  that 


(5.4) 


(5.5) 


li<I,*)|  <  ud  +  IH1+C) 

+  lll5> 


If  Y  and  0  are  small  enough,  then  (5.2)  has  a  unique  solution  of  the  form 


(5.6) 


h(z)  -  3<j>(z)/3g  +  hUH*) 


where  L  (B)  a  A  ►  k(A)(*)  e  L_(8)  is  continuous  and  compact. 

—  a  — •  a  . 

Proof.  It  suffices  to  take  8  *  B^.  We  use  the  following  amalgamation  of 
results  of  Sobolev  and  Kantorovich  (cf.  Sobolev  1950,  5®)  *nd  Kantorovich  & 
Akilov  (1977,  Chap.  XI,  §3):  If.  f  €  L« (B  ) ,  with  S  >  1,  then  there  is  a 


S  >vy*, 

ST  --  % i 

\*S  -  * 
V  -J 


VV>j 

n 

■  '  •  '(L.  1 

-  -V  V 
w*  /  r. 

( v 

s  *»'  v 

V*. 
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» -V  .V 
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continuous  function  Y  ♦  ie(8»Y)  that  strictly  increases  from  0  to  •  as_  Y 
increases  from  0  to  •  such  that 

(5.7a)  lKf,Lv(BY)l  <  K(6,Y)lf,Lg(BY)l 

where 

f(x) 

(5.7b)  (Kf ) (z)  =  /  - r  dv(z)  , 

B  |x  -  z|2 

Y  1  IS*  fW1 

v  -  •  i B  >  3,  v  <  3B/(3  -  8)  if.  8  <  3.  Moreover  K  is  compact  (and 

continuous)  from  L„(B_.)  to  L.(B..). 

————————  p  Y  —  v  y 

We  wish  to  show  that  &  *  &(h,£» * )  is  a  contraction  from  L^B^)  to 

itself.  Let  h  £  Lq(By).  We  first  identify  f  of  (5.7)  with  (the  components 
*  -1 

of)  and  chose  8  ■  o(1  +  C)  •  (Then  8  <  3  if  and  only  if 

A  -1 

a  -  35  <  3.)  Then  (5.4)  ensures  that  1  £  L  (B  ).  Since  38(3  -  8)  >  a  when 

'*•  oi  Y 

8  <  3,  we  can  take  v  »  u  in  this  case.  Thus  (5.7)  implies  that 
h  ►  kth^)  (•)  maps  L^B^)  into  itself. 

We  now  show  that  fc  *►  2s(fe#£,  •)  is  a  contraction#  Let  c  **a*BY*' 

6fe  5  fe1  -  fe2*  Then 

(5.8)  l«)s(brb2.?H*>!  s  |Js(fcrI>(s>  -  Js(fe2'I)(^l 


+  (1 


-  t)h2(x),£(x),x)| |6h(x)j 


dv(x) 


where  t  c  10,1].  We  now  identify  f (x)  with  the  numerator  of  the  integrand  in 
the  right-most  terra  of  (5.8).  We  henceforth  suppress  the  arguments  of  the 
functions  appearing  in  this  numerator.  Let  us  choose  8  and  v  as  above. 


noting  that  8  <  a.  From  (5.7),  (5.8),  and  the  Hfilder  inequality  we  then  obtain 


(5.9) 


^JS'W*  <  •c<b*y>C/  |ai/ah|8|«h|Bdv>vs 

N 

<  K(e,Y)»3y3h,LaB/(a_8)(BY)ll«h,La(BY)l  . 

A 

Since  o$/(a  -  0)  -  B,  condition  (5.5)  implies  that  3}/3h  €  ^oB/fo-B ) ^y ^  80 

that  the  rightmost  term  of  (5.9)  is  well  defined. 

A 

We  now  prove  the  compactness  of  ^  »  £(&,#*).  The  compactness  of  K 
introduced  in  (5.7b)  implies  that  the  mapping 

r  ,  *  (*  -  *>  , 

(5.10)  - 3 - dv(s)J  e  La(BY) 

A 

is  compact.  Condition  (5.4)  ensures  that  (h,£_)  *  <J(h»5.»*)  takes  La(By)  to 
Lg ( 3^ ) .  By  the  properties  of  Nemytskii  operators  (cf.  Krasnosel ' skii  (1956, 

Sec.  1.2) ,  this  mapping  is  continuous.  Since 

A  A 

,  4<H(*)/3j5  +  JS(5,X),E.(X))S)  -  (»  -  *)«»v(x)  - 

La(B  )3AH.J - ; - -3 - 5 

®Y  Is  ~  JS I 

is  the  composition  of  a  compact  with  continuous  operators,  it  is  compact.  o 
It  follows  from  (5.5)  and  the  properties  of  k  that  we  can  make  the 
coefficient  of  !6h,Lal  in  the  right  most  term  of  (5.9)  less  than  unity  by 
fixing  0  and  taking  y  small  enough  or  by  fixing  y  and  taking  0  small 
enough  or  by  taking  each  small  enough.  The  Contraction  Mapping  Principle  ensures 
that  (5.2)  has  a  unique  solution  giving  h  as  a  continuous  function  of  A_.  The 
composite  function  obtained  by  substituting  this  solution  h  into  £(£,£,*)  has 

A 

value  denoted  by  ^{A^x) . 

A  number  of  related  results,  including  some  for  unbounded  domains,  can  be 
based  on  the  techniques  presented  by  Sobolev  (1950,  §56,9),  Stein  (1970,  Ch.  V), 
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and  Kantorovich  &  Akilov  (1977,  Chap.  XI).  Note  that  Theorem  5.4  says  that  (5.6) 

A 

is  valid  provided  the  dependence  of  ^  on  h  becomes  weaker  as  the  body  B 
becomes  larger. 

A 

He  now  suppose  that  ^  is  such  that  admits  a  representation  of  the  form 
(5.6).  He  substitute  (5.6)  into  the  right  side  of  (4.2)  to  get 
(5.11)  T(z)  -  T(F(s),g(r),A(r),3(p(E)/az,3*(z)/3r  +  k(4,j),x)/  etc. 

Our  governing  equations  are  obtained  by  substituting  (5.11)  into  (3.5)- 
( 3.8 ) t 


A 


(5.12) 

Div  T  -  0  , 

(5.13) 

Div  9+3* (84/8*) 

(5.14) 

A 

Div  d  -  0  , 

A 

(5.15) 

Div  b  -  0  , 

A#  r 

where  the  arguments  of  the  constitutive  functions,  decorated  with  carets,  are 
indicated  in  (5.11).  Equations  (5.12)-(5.15),  having  six  scalar  components ,  form 
a  quasilinear  system  of  partial  functional  differential  equations  for  the  six 
unknown  components  of  y,A,(p,\|i.  All  other  variables  we  have  introduced  can  be 
expressed  in  terms  of  these. 

#w  A# 

Toupin  (1956)  took  the  polarization  g  and  magnetization  m  as  independent 
constitutive  variables.  One  of  the  goals  of  our  paper  is  to  exhibit  the 
mathematical  advantages  of  choosing  d,b , ^  to  be  dependent  constitutive 
variables  and  choosing  <§  and  Jj  to  be  independent  constitutive  variables.  In 
this  regard  we  generalize  formulations  of  Pao  &  Hutter  (1978),  Jordan  &  Er ingen 


(1964),  and  Ersoy  &  Kiral  (1978)  inter  alia 


6.  Ellipticity  and  Growth  Conditions 

Our  basic  constitutivs  assumptions  are  expressed  in  terms  of  the  quadratic 


(6.1)  2 

V(3i/3£>tft  +  &*OS/3g)*fc  +  &«OT/3e)*a  +  A«Os/3h)*v 

+  £*<3a/a£>!A  +  Og/3jj)  *t  +  t» (3g/3e)  *u  +  t«(3g/3h)«v 

A  A  A  A 

+  u»(3d/3F):A  +  u*(3d/3jj)*t  +  u»(3<J/3e)»u  +  u»(3d/3h)»v 

A  A  A  A 

+  v»(31j/3f)iA  +  v*  (3Jj/3jj)  *t  +  X*(3b/3e)»u  +  v*  ( Sb/3h)  * v  . 

A  A  A  A 

If  >  0  V<A*;6»H»X>  ^  then  (T,a,d,b)  i8  8aid  to  be 

strictly  monotone.  The  use  of  this  attractive  mathematical  restriction  would 
deprive  the  theory  of  much  of  its  physical  versatility.  Among  its  adverse 
consequences  (discussed  in  detail  by  Antman  (1983))  is  that  the  uniqueness 
theorems  it  implies  effectively  prevent  the  buckling  of  a  column  of  such  a 
material  however  slender  under  a  compressive  load  however  large. 

We  can  eliminate  this  kind  of  uniqueness  in  the  mechanical  response  while 
preserving  it  fully  in  electromagnetic  response  and  partially  in  the  thermal 
response  by  weakening  the  strict  raonotonicity  condition.  If  o>(rs,t,u,v)  >  0 

A  A  A  A 

then  we  8*y  that  (T.g.d.b)  satisfies  the  (strict 
form  of  the)  restricted  strong  ellipticity  condition.  (Note  that  IV  equals  the 
dyadic  product  ££  if  and  only  if  ft  has  rank  1 . )  In  Section  8b  we  discuss  the 
physical  significance  of  the  restricted  strong  ellipticity  condition.  If 

A  A  A  A 

>  0  v(t***5l'*2'53)  **  then  (T,q,d,b)  satisfies 

the  (strict  form  of  the)  strong  ellipticity  condition.  This  condition  is  the 
generalization  to  elliptic  systems  in  divergence  form  of  the  Legendre-Hadamard 
condition  of  the  calculus  of  variations. 

In  this  paper  we  shall  study  the  strong  ellipticity  condition  and  its 
restricted  form.  Our  subject  is  insufficiently  developed  to  determine  whether 


w-vVV 


phenomena  permitted  by  the  strong  ellipticity  condition,  but  prohibited  by  its 
restricted  form,  are  observed  (cf.  Sec.  8b).  The  intuitive  content  of  these 


conditions  is  described  in  Section  1.  Most  special  theories  of  material  behavior 
of  electro-magneto-thermo-elasticity  satisfy  the  restricted  strong  ellipticity 
condition  because  many  of  the  "off-diagonal"  terms  in  (6.1)  are  zero.  (But 
recent  work  on  the  study  of  plastic  effects  and  phase  changes  treates  theories  of 

A 

elastic  solids  for  which  even  rs: ( 3T/3F) srs  need  not  be  positive  for  jjg  g. 
Cf.  Ericksen  (1980).) 

we  now  study  the  behavior  of  the  constitutive  equations  at  extreme  values  of 
their  arguments.  The  conditions  we  impose  must  be  consistent  with  the  strong 
ellipticity  condition.  Since  our  work  is  just  a  generalization  of  that  of  Antman 
(1983),  we  omit  an  extensive  commentary.  In  Sections  8  and  11  we  describe  more 
specific  conditions  appropriate  for  special  problems. 

Recall  that 

(6.2)  i  s  x • 

Let  £  and  c  be  unit  vector  fields  depending  on  X'X'S*  strong 

ellipticity  condition  implies  that 


(6*4)  P(&s)  -  {j&*£*S  «  Rjdet  £  >  0} 

is  eithsr  an  opan  half-line,  or  the  whole  line,  or  empty  and  can  then  be  written 

as 


(6.5) 


P(ag)  -  (A  (ft£)#*  (*c))  . 


We  suppress  the  dependence  of  V  and  A  on  £,£,*.  The  facts  motivate  the 
following 

6.6.  Hypothesis.  Let  a  and  c  be  unit  vector  fields  depending  on  £,;£,*.  .If 

A 

tXfig)  is  an  open  half-line  or  the  whole  line  and  if  a*F»c  *  a»T»c  is  strictly 
increasing,  then 

A  <«  X 

(6.7a)  a*T*c  ♦  +•  as  a*F»c  +  A  (ac)  for  fixed  T  -  ((a»F»c)ac,0,0,0,0),v,z  . 
If  jj*a  w  <j»a  is  strictly  increasing,  then 

(6.7b)  3* a  ♦  ±*  as  £•£  ♦  +"  for  fixed  £  -  (0,(g*a)a,0,0,0),y,z  . 

A 

If  e»a  ►  d»a  is  strictly  increasing,  then 

A 

(6.7c)  ♦  ±-  as  ♦  +“  for  fixed  J_  -  (g#0#0 ,  (e»a)a,0 )  ,y,s  . 

A 

If  h*a  *■  b*a  is  strictly  increasing,  then 

A 

(6.7d)  ♦  +•  as  ♦  ±“  for  fixed  £-  (0,0, 0,0, (h*a)a) ,y,z  . 

If  1X&£)  is  a  half-line,  which  happens  exactly  when  the  cofactor  of  a*£*c 
in  det  F  does  not  vanish,  then  3P(jjg)  is  a  point  (either  A+(ac)  or 
A  (ac)).  Then  (6.7a)  implies  that  we  can  define  a  function 

(6.8)  eg  ►  6(ac)  e  {— 1 , 1 > 
such  that 

A 

(6.9)  6(ac)a«T»c  ♦  -•  as  a*F*c  ♦  9P(ac) 

MM  M  M  M  M  M  M  MM 

for  fixed  £  -  ( (a*F*c)gc,0,0,0,0) ,y,z  .  Note  that  det  F  ♦  0  so  the  local 
volume  ratio  shrinks  to  0  as  a»F*c  ♦  3P(ac). 

M  M  M  MM 

We  now  complement  Hypothesis  6.6  in  a  way  that  promotes  the  analysis  of  Part 

A 

II  by  describing  the  behavior  of  the  constitutive  function  T  as  more  than  one 
component  of  £  are  allowed  to  vary. 


m, 
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,...,9}  be  a  basis  for 


Lin  consisting  of  dyadic  products  of  unit  vectors,  which  may  depend  on  r.,x,jg. 
and  let  {et>  be  the  basis  dual  to  {jgT}.  Let  F:^  be  strictly 

increasing  for  each  t .  Let  {J^}  be  a  sequence  of  states  such  that  the 
l  (E^ )  formed  from  {r^}  are  actually  independent  of  n.  Let  the  set  of 


integers 

{l,...,9}  be  written  as  a  disjoint  union 

a  U  b 

ucudueu^  with 

(6.11a) 

*  * 

3P(E)  t  <>  and  6<ET)T<r  ,z)  :E  —  for 

T  C  CL 

9 

(6.11b) 

3P(EX)  f  *  and  £n:;gT  +  3I>(gT)  for  r  e  b 

t 

(6.11c) 

3P(E_)  ?  6  and  F  :E_e  compact  subset  of 

)#A+(E  ))  for  tec. 

(6. lid) 

3P(E)  +  for  red, 

- 

l 

T 

(6. lie) 

3P(E_)  =  6  and  F  :E_e  compact  subset  of 

),i+(E  ))  for  T  €  e  , 

(6. Ilf) 

3P(gT)  =  6  for  t  e  . 

T  — — 

If  a  + 

,  then 

(6.12a) 

6(E_ )T( T  ,z):E_  +  -“  Vx  e  a  U  b  U  Cl  for  each  t 

in  d  either 

*  • 

If  :E  ♦  •  or  6(E_)T(T  ,z):E,  +  -®. 

If  b  f1  fb,  then 

A  * 

(6.12b)  either  (i)  6(E,)T(T  ,z) :E_  ♦  -«  VT  e  d  U  d  U  c  U  d  or  else 

■  1  ■  ■  -  t  ~  — n  ~  ~x  ■  ■  ■ 

(ii)  3t  e  CL  U  d  such  that  |f  :E_|  ♦  »  and  3t  ta  Ufc  Uc  Urf 

-  i~n  '  - 

such  that  |t(T  ,z):E  |  ♦ 

Moreover,  the  dualization  obtained  by  respectively  replacing  the  statements 

A  * 

F  :E„  ♦  3P,(E„),|F  :E_  I  +  •,  fiT(r  ,z):E,  ♦  — 


appearing  in  (6.11),  (6.12)  by  their  opposites 


F  :E  ♦  3P (E  ) ,  5T( T  ,z)sEt  +  • 
~n  ~t  — T  ~  -n  ~  ~x 


is  also  valid. 

The  statement  containing  (6.12b)  may  be  loosely  interpreted  thus:  If  there 
are  fibers  compressed  to  zero  length  in  some  directions,  then  either  the  material 
is  squeezed  out  with  an  infinite  stretch  in  another  direction  or  else  it  is 


-20- 


prevented  from  doing  so  by  infinite  compressive  stresses  in  all  other  directions. 
The  other  statements  have  similar  interpretations.  The  whole  hypothesis 
effectively  says  that  extreme  behavior  in  one  direction  must  be  accompanied  by 
extreme  behavior  in  some  other  direction.  This  observation  is  used  in  Part  II  to 
establish  regularity  results  by  showing  that  behavior  could  be  extreme  in  only 
one  direction  and  therefore  cannot  be  extreme.  It  would  be  easy  to  generalize 
Hypothesis  6.7  to  account  for  extreme  couplings  between  the  mechanical, 
electromagnetic,  and  thermal  effects,  but  the  intuitive  evidence  for  such  a 
generalization  is  not  compelling. 

It  is  important  to  note  that  the  transformations  (3.3)  and  (3.4)  ensure  that 
the  Maxwell  stress  contains  terms  with  (det  E)'1  as  factors.  These  terms  could 
compete  with  the  "purely  mechanical  part  of  this  stress"  when  det  E  is  small. 
Our  constitutive  hypotheses  on  the  effective  stress  control  this  competition. 

They  say  that  the  material  response  in  large  compression  is  dominated  by  that  for 
purely  mechanical  response. 


7.  Boundary  Conditions.  The  Principle  of  Virtual  Work 


At  a  boundary  point  z  €  38  we  may  prescribe  the  position  y(g)  or  merely 
subject  it  to  certain  constraints,  such  as  the  requirement  that  it  be  confined  to 
a  fixed  surface.  To  account  for  the  varied  possibilities  it  is  convenient  to 
describe  such  boundary  conditions  in  the  language  of  holonomic  constraints.  We 
accordingly  specify 

(7.1)  %(z)  «■  2(z,r)  for  each  z  €  3B 

where  X  is  »  given  function  continuously  differentiable  in  _r  e  R3,  which 

represents  the  set  of  generalized  coordinates  for  y(g).  ®ie  rank  of  ^X^IL  i® 

the  number  of  degrees  of  freedom  of  g.  Equation  (7.1)  restricts  ytg)  to  a 

#  -  # 

manifold.  The  set  of  vectors  X  ^  of  form  *£.  for 

r*(z)  €  R3  form  the  tangent  space  to  this  manifold  at  y(z).  The  elements 
X  (j)  of  this  tangent  space  are  called  virtual  displacements.  We  complement 

(7.1)  by  specifying  the  projection  of  the  traction  T«n  on  this  tangent  space: 

(7.2)  *2*2*  "  t)*  (Sjg/Sr)  ■  JO  at  each  z  e  38 

where  t  is  a  given  function  of  y'g),  A(g),  »(g),  ♦  (&)*  z  and  possibly  other 
variables.  Thus 

(7.3)  (n*T*  -  t)*y#  =0  on  3B  . 

<V  A#  r>j  Aj 

At  each  2  c  38  we  also  prescribe 

(7.4a,b)  either  A(z)  *  K(z)  or  g(z)*n(z)  =  Y(£(z)  *  A(z)  *<l>(z) *♦{*>  »5>  » 

(  7.5a,b)  either  v(z)  *  ?(z)  or  d(z)»n(z)  =*  S(;g(z) ,  A(z)  ,<p(z)  ,i|>(z)  ,z)  , 

(7.6a,b)  either  <l<(z)  =  iji(z)  or  b(z)*n(z)  *  ,\[z)  ,<r(z)  ,  z)  . 

Let  A  ,ij»  be  arbitrary  continuous  functions  on  3B  that  respectively  vanish 
where  (7.4a),  (7.5a),  (7.6a)  hold.  They  are  virtual  fields.  Then  in  analogy 
with  (7.3)  we  have 

(7.7)  (g»n  -  y)A*  +  (d*n  -  J)<p*  +  ( b*n  -  8)4>*  =  0  on  38  • 


Our  fundamental  equations  of  balance  are  the  integral  versions  of  (3.5)- 
(3.7),  which  are  to  hold  over  "almost  all"  nice  subbodies  of  B.  These  equations 


can  be  supplemented  by  appropriately  weakened  forms  of  the  boundary  conditions  we 
have  just  listed.  Antman  a  Osborn  (1979)  show  (strictly  speaking,  for  the  purely 
mechanical  problem)  that  when  all  the  integrals  make  sense  as  Lebesgue  integrals, 
then  these  equations  and  boundary  conditions  are  equivalent  to  the  Principle  of 
Virtual  Work 

(7.8)  /  [T:Ox#/3a>  +  g*(3A#/3*)  +  d»(3*#/3a)  +  b- ( 3<|>#/3 z) ] dv 

8 

r  #  *_#  „  # 

-  J  (f*x  +  PA  -  o<p  )dv 

8 

“  /  +  +  8<»#)da 

38 

for  all  reasonably  nice  fields  x  '7  *♦  having  the  boundary  behavior 
specified  above.  Equation  (7.8)  is  just  the  weak  formulation  of  our  boundary 
value  problem  consisting  of  ( 5. 12)-( 5. 15)  subject  to  (7.1),  (7.2),  (7.4)-(7.6). 

0k 

The  arguments  of  T ,  etc . ,  are  given  in  (5.11). 

In  many  circumstances  the  deformation  of  a  body  subjected  to  the  action  of 
external  electromagnetic  fields  changes  the  ambient  fields.  Thus  there  would  be 
a  complete  coupling  between  the  fields  interior  and  exterior  to  the  body.  Since 
our  goal  is  to  study  the  role  of  the  constitutive  assumptions  of  Section  6,  we 
are  avoiding  such  coupled  problem  by  restricting  electromagnetic  boundary 
conditions  to  (8.5)  and  (8.6).  Methods  for  treating  fully  coupled  problems  would 
be  similar  to  those  of  Section  5. 


8.  General  Existence  Theorems 

In  this  section  we  obtain  existence  theorems  for  two  important  special 
classes  of  problems ,  which  can  be  readily  treated  by  means  of  recent  results  for 
elliptic  systems.  For  the  first  problem  we  assume  that  there  is  neither  thermal 
nor  electrical  conduction  and  that  there  is  a  stored  energy  function}  thus  this 
reduced  problem  admits  a  variational  formulation.  For  the  second  problem  we 
assume  that  the  material  is  rigid.  The  restricted  strong  ellipticity  condition 
then  reduces  to  a  monotonicity  condition,  which  is  capable  of  handling  our 
nonlocal  operators. 

a.  Conservative  Problems 

We  assume  that  the  material  does  not  conduct  electricity  so  that  the 
constitutive  function  jJ,  =*  0.  Thus  the  Joule  heating  is  zero  (cf.  (5.13)). 
Moreover,  (5.2)  reduces  to 

(8.D  &  -  . 

We  assume  that  one  of  the  following  conditions  holds: 

i)  3*P»Y  depend  only  on  g,X,z  and  the  boundary  value  problem  (5.13),  (7.4) 
has  a  weak  solution  X  in  a  suitable  Sobolev  space.  (In  part  (b)  below,  we  show 
how  a  slight  strengthening  of  our  hypotheses  ensures  the  existence  of  X.) 

ii)  The  boundary  value  problem  (5.13),  (7.4)  has  a  solution  X  independent  of 
the  fields  F,  3?/3z,  3i|>/3z.  This  situation  would  occur  if  3*0  when  g  *  0 
and  if  X  is  prescribed  to  be  constant  X^  on  38,  for  then  the  boundary  value 
problem  would  admit  a  solution  X  =  X^  on  all  of  8. 
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iii)  The  constitutive  functions  T,d,b  are  independent  of  g  and  X. 

In  cases  (i)  and  (ii)  we  can  substitute  the  solution  X  and  its  gradient 
g  into  (5.12),  (5.14),  (5.15),  (7.1),  (7.2),  (7.5),  (7.6).  Since  X  is  known. 


its  presence  in  these  equations  merely  changes  the  dependence  of  the  constitutive 
functions  on  £.  Zn  case  ( iii) ,  these  equations  are  unaffected  by  the  solution 
X. 

We  assume  that  c  is  a  prescribed  function  of  j.  We  finally  assume  that 
there  is  a  stored  energy  function  W  depending  on  £*g*X,g,^,jg  with  W 
continuously  differentiable  in  £,g,b,  continuous  in  g  and  X,  and  measurable 
in  £,  for  all  values  of  the  remaining  arguments ,  such  that 

(8.2)  £  -  3W/3F,  d  -  3W/3e,  b  -  3w/3h  . 

(The  Clausius-Duhera  inequality  would  deliver  a  specific  thermodynamic  function 
for  W  and  show  that  it  would  be  independent  of  g.)  The  discussion  following 
assumptions  (i) ,  (ii),  (iii)  motivates  us  to  suppress  the  dependence  of  w  on 
g  and  X,  their  effects  being  absorbed  by  the  dependence  of  W  on  g. 

We  suppose  that  the  body  force  £  is  conservative  so  that  there  is  a 
function  (y,z)  U(y,j6),  with  0(*#jj)  continuously  differenciable  for  all  g 

in  B  and  with  U(y,  •)  measurable  for  all  y  e  S3,  such  that 

(8.3)  l  ■  ”30/3y  . 

We  suppose  that  t,B  ,B  of  (7.2),  (7.5),  (7.6)  are  conservative  so  that  there  is 
a  function 

I3  *  I3  x  |3  *  38  j  (y ,e,h,z)  w  V(y,e,h,z) 

a#  a#  a#  n#  a.  ax  ax 

with  V( *,*,*,g)  continuously  differentiable  for  all  £  in  38  and  with 
V(y,g,^,«)  measurable  for  all  y,£,tj,  such  that 

(8.4)  £  -  3v/3y,  J  -  3v/3g,  5  -  3V/3h  . 

(The  domain  of  V(y,g,ft,*)  may  be  taken  to  be  the  closure  of 

38\{g  €  38*y*(z)  •  0,  e*(z)  "  0,  <{i*(z)  ■  0}.  (See  Section  7.)  V  could 
conceivably  depend  on  g  and  X.  We  suppress  any  such  dependence  in  accord  with 
the  policy  we  have  adopted  above. 


We  finally  assume  that  38  is  bounded  and  has  a  locally  Lipschitz 

#  # 

continuous  graph.  Moreover,  we  require  that  the  supports  of  <p  ,i|»  ,  and  the 

* 

components  of  x  be  nice  enough  to  ensure  that  the  boundary  conditions  (7.1), 
(7.5a),  (7.6a)  are  assumed  in  the  sense  of  trace  when  lie  in  Sobolev 

spaces  of  the  form  W^(B)  with  p  >  1.  (Necessary  conditions  for  these 
properties  are  not  known.  See  the  discussion  of  Antman  &  Osborn  (1979).) 

Under  these  conditions  the  weak  form  of  the  Euler-La grange  equations  for  the 
functional 

(8.5)  X(x»'f+)  “  I  (*>#  f|  (*>#  §£  (*>»*)  +  U(X(*),z)  +  0(E)*(z)]dv(*> 

B  ~  ~  ~ 

+  /  V(x(2>/*(2>/*(2>/*><Ja(z) 

3B 

for  satisfying  (7.1),  (7.5a),  (7.6a)  have  exactly  the  form  of  (7.8)  with 

# 

X  =  0.  (Of  course,  many  authors  take  a  variational  principle,  such  as  this,  as 
the  starting  point  for  the  derivation  of  the  governing  equation  for 
electromechanical  interactions.  See,  e.g.,  Toupin  (1956),  Brown  (1966),  Nelson 
(1979).) 

x 

Let  F  denote  the  cofactor  tensor  of  J.  W  is  said  to  be  polyconvex  (cf. 
Ball  (1977))  if  it  can  be  written  in  the  form 

W(F,e,h,z)  -  (l(F,F*,det  F,e,h,z) 

a/  a/'  »w  (v  w  w  n 

3  9 

with  fl( •,•»•,•  ,*»;g)  convex  on  Lin  *  Lin  x  (0,*)  x  s  x  tr  for  each  z  e  8. 

The  work  of  Ball  (1977)  shows  that  if  W  is  polyconvex,  then  (10.2)  satisfies 
the  restricted  strong  ellipticity  condition  of  Section  6.  To  account  for  (6.7) 
we  require  that 

(8.6)  Q(F,F*,5,e,h,z)  ♦  »  as  6  0  . 

8.7.  Theorem.  Let  W  be  polyconvex  and  satisfy  (8.6).  Let  there  be  numbers 


X  e  L^(38)  such  that 

(8.8)  >  “(g)  +  *(IeIp  +  lgxlq  +  «r  +  !s!p  +  lfelp) 

for  all  jg  €  8, 

(8.9)  tf(X*£)  >  “(g)  *or  a11  J6  «  8  , 

(8.10)  a  €  ^/(p_i)(8)  i 

(8.11)  VlX/P^/g)  >  X(*)  for  all  3  e  3B  . 

Let  38  have  the  properties  specified  above.  Let 

(8.12)  W  2  {(x,»,*>  €  Wp(8jJX  C  Lq(8>,  det  F  €  I^fB)  , 

det  F  >  0  a.e.  for  F  -  &x/3*»  (7.1),  (7.5a),  (7.6a) 

are  satisfied  In  the  sense  of  trace  where  they  are  prescribed 

on  3B  • } 

If  there  exists  an  element  (x1»9^»’l,1)  «  W  such  that  I(j^»9^,^)  <  •,  then 
there  exists  an  e lenient  (y ,?,3»)  that  minimizes  I  on.  W. 

The  proof  of  this  theorem  is  effected  by  making  minor  adjustments  to  that  of 
Ball  (1977)  and  is  accordingly  omitted.  (Further  developments  of  Ball's  theory, 
useful  for  our  class  of  problems,  are  given  by  Ball  &  Murat  (1984),  Dacorogna 
(1982),  and  the  references  cited  therein.) 


b.  Rigid  Conductors 

We  now  study  the  effects  of  the  conduction  of  heat  and  electricity,  but 
confine  our  attention  to  rigid  bodies,  for  which  F  is  constrained  to  be  the 
identity  I.  We  accordingly  take  the  virtual  displacement  x  »  appearing  in 
(7.8),  to  be  0.  This  choice  ensures  that  the  First  Piola-Kirchhoff  stress 
tensor,  which  now  is  the  Lagrange  multiplier  maintaining  the  constraint  of 
rigidity,  does  not  enter  into  (7.8).  We  drop  F  from  the  list  of  variables 
constituting  in  (5.1)  and  from  the  arguments  of  the  constitutive  functions 


S'3'fe'i  (c£.  (5.11)).  Our  boundary  value  problem  reduces  to  (5.13)-(5. 15), 
(7.4)— (7.6),  whose  weak  form  is  the  suitably  specialized  version  of  (7.8). 

We  assume  that  the  restricted  strong  ellipticity  condition  holds.  Thus 
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<jg*S*h)  **■  >  is  strictly  monotone. 

This  condition  prohibits  certain  kinds  of  nonuniqueness. 

Since  hysteresis  frequently  is  associated  with  nonuniqueness  and  since 
hysteresis  is  one  of  the  most  important  phenomena  of  ferromagnetism,  it  might 
appear  that  our  use  of  the  restricted  strong  ellipticity  condition  precludes  us 
from  dealing  with  ferromagnetic  materials.  But  molecular  theories  of 
ferromagnetism  (cf.  Tebble  (1969))  suggest  that  hysteresis  is  associated  with 
constitutive  equations  with  nonlocal  effects.  If  we  accept  such  theories,  then 
to  account  for  ferromagnetism  it  is  necessary  to  generalize  the  form  of  our 
constitutive  functions  before  relaxing  the  ellipticity  conditions.  We  do  not 
attempt  such  a  generalization  here:  Our  analysis  should  be  regarded  as  merely 
applying  to  paramagnetic  materials.  We  do,  however,  examine  nonlocal  operators 
that  are  introduced  by  the  mathematical  approach  we  use  to  handle  electric 
currents.  Some  of  the  methods  we  use  can  be  applied  to  more  general  kinds  of 
nonlocal  behavior. 

We  now  outline  an  existence  theory  that  can  be  applied  directly  to  our 
specialized  version  of  (7.8).  We  first  present  the  theory  in  an  abstract  form  in 
order  to  facilitate  a  comparison  of  it  with  presentations  in  the  mathematical 
literature.  Afterward  we  make  the  requisite  identifications. 

Let  8,  as  before,  be  the  closure  of  a  domain  in  R  .  We  assume  that  SB 
has  a  locally  Lipschitz  continuous  graph.  A  typical  point  in  B  is  denoted  z. 
Let  u(z)  *  (u„ (z) , . . . ,u  (z)).  For  p  €  (1,w),  let  the  operator 

—  I  ~  HI  ^ 

(8.13)  L  (B)m  x  L  (B)3m  3  (u, 3v/dz)  w-  k(u,3v/3z> ( • )  e  [L_(8)]r 

P  P  —  —  <v  —  —  —  P 

take  bounded  sets  into  bounded  sets.  Let 

3m  r  €  R 

(8.14)  *  R  *  r  *8  3  *  -,  i  -  1,...,m  , 

B1(C,n,£,z)  €  R 

(8.15)  if  *  38  3  (£,z)  »  y1(C.,z)  e  R,  i  -  1,...,ra 

satisfy 
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(8.16a)  For  almost  all  jj  in  8,  the  functions  aS •,»,•,£),  B^t •,•,•,*)  are 
continuous  and  for  all  £,.n.,£>  the  functions  a*1  (£,£,£,  • ) ,  B*”(£,n,£,») 
are  measurable.  (These  are  the  Carathftodory  conditions.) 

(8.16b)  For  almost  all  g  £  38,  the  functions  yN*,*)  are  continuous  and  for 
all  £,  the  functions  Y1^,  •)  are  measurable  on  3B  (with  respect 
to  two-dimensional  Lebesgue  measure). 

(8.16c)  There  exist  a  constant  c-j  >  0  and  a  function  k^  e  Lp*(8)  (with 
p*  ■  p/(p  -  1))  such  that 

|ai(C.,n_,£,*)  | ,  |6i(C,n,£,z)|  <  cn  C  |  C  |  1  +  M*”1  +  Id*"1  +  k^z)] 

for  i"  1,...,mf  o  ■  1,2,3. 

The  Holder  inequality  then  implies  that  the  functions 


i  3v  *  3u  .  3v  »  3u 

ft  <«<*>'  ^  f>.  8  (u,  (•),  jt(u,  ajH*),*) 


are  in  I*p*(8)  for  all  u,;/  £  (Wp)™.  It  follows  that  the  functional 

j  3u  «  3u  3w 

(8.17)  a(u,w)  5  /  tl  ft  <£<£>'  -§7  <ft>*  IjlH.*  a£> <£>'£>  *  37*  <*> 

8  i  ~  ~  ~ 

3u  «  3u 

+  l  BNu(z),  -jy  (z),  k(u,  ■5Y)(z),z)wi(z)]dv(z) 


+  /  I  Yi(u(z) )w. (z)da(z) 
38  i  “  1 


is  well  defined  for  all  u.,w.>  e  Wp(B). 

We  shall  prescribe  u^x.fUg,  respectively  on  subsets  Sf,...,Sm  of  38. 
We  assume  that  these  subsets  are  measurable.  Let  V  be  the  closed  subspace  of 
[Wp(8)]ra  containing  (W^(8)]°  that  consists  of  functions  (w1,...,wm)  for 
which  w^  ■  0  on  S1,...,wm  *0  on  Sm  in  the  sense  of  trace.  Let  u  be  a 


given  element  of  [Wp(8)]  .  We  require  that  agree  with  on  S  j ,  etc., 

in  the  sense  of  trace  by  seeking  solutions  u.  of  our  equations  in  [Wp(B)]m  to 
which  \i  -  u  €  V.  (This  prescription  of  boundary  conditions  enables  us  to  avoid 
the  very  delicate  questions  of  whether  functions  defined  on  S^,...,Sm  can  be 
extended  to  functions  in  Wp ( 8 ) . ) 

Since  V  a  w_  ►  a(_u»w_)  is  a  bounded  linear  functional  for  each  in 
Wp(B),  the  Riesz  Representation  Theorem  enables  us  to  write 

(8.18)  a(u^w)  *  <M_u) ,  w>  , 

where  <_v,w>  =  /  v*w  dv.  If  the  a1  are  continuously  differentiable  and  if 
8 

is  twice  continuously  differentiable  on  8  and  vanishes  on  3B,  then 
A(u)  -  A1  (u) , . . . ,  Ara(_u) )  where 

3u  A  3u  3u  3u 

(8.19)  aMu)  =  -Div  aNu,  k(u,  -j^,z)(z),z)  +  B(u#  ^,lt(u,  -j^-,z)  (z)  ,z)  . 


Let  us  set  *  (J^,..  0ur  basic  abstract  result  is  the  following: 

8.20.  Theorem.  Let  38  have  a  locally  Lipschitz  continuous  graph.  Let 
p  €  (1,®).  Let  (8.16)  hold.  Suppose  that 

a(v,v) 

(8.21)  ~Tv~Vl  *  °°  -2£  00  for  v,  c  V  , 

(8.22)  I  ,rj.  +  P,£,z>  -  >  o  VP  5  (j^,...,^)  f  0_  . 

i 

If  8  is  bounded,  let 

-1 

(8.23a)  l  ai(C.,n,C,z),n.  tlni  +  |nJP  J  ®  as  |nj  -►  ® 

i  ~  ““  1 

for  almost  all  g  in^  8  and  for  bounded  lf_  8  is  unbounded,  let  the 


following  stronger  restriction  hold!  There  Is  a  number  C2  >  0  and  a  function 
k2  €  1^(8)  such  that 

(8.24b)  l  *t(£#2!>#£#J8>,3i  >  e2|n|p  -  k2(z)  . 


Define  k  by 

(8.25)  [Wp(B)]B  a  u  **•  k(u)(»)  =  k(u,3u/3z)(*) 

A 

where  It  is  defined  in  (10.13).  Let  be  the  characteristic  function  of  a 

set  C  ^n  8.  For  every  subdomain  C  of  8  with  compact  closure  in  B  let 

(8.26)  [W^(B)]“  9  u,  *►  Xc(*)k(u)(* )  €  [Lp(C)]r 
be  compact.  Then  for  every  _f_  £  V*  and  for  every  u  €  [Wp(8)lm  there  exists 
a_  _u  e  [Wp(8)]m  with  u_  -  u  e  V  such  that 

(8.27)  <A(u)  ,v>  -  <t_,v>  Vv  €  V  . 

The  proof  of  this  theorem  is  obtained  by  making  minor  adjustments  to  those 
of  Brezis  (1968)  (cf.  Lions  (1969,  p.  297))  and  Browder  (1977).  We  note  the 
following  points t  In  a  bounded  domain  the  operator  A  is  of  the  "calculus  of 
variations  type"  because  of  its  monotonicity  in  the  local  values  of  its  highest 
order  derivatives  and  because  of  its  compactness  (through  k)  in  the  global 
values  of  the  highest  order  derivatives.  Since  our  integral  operator  (7.6)  for 
constant  electric  currents  is  not  compact  on  unbounded  domains,  we  had  to  use  the 
theory  of  Browder  (1977)  based  upon  the  compactness  of  (8.26)  to  support  our 
intended  applications. 

We  identify  the  variables  appearing  in  Theorem  8.20  with  those  used  in  the 
problem  outlined  at  the  beginning  of  this  subsection.  In  particular,  we  set 

(8.28)  u,  -  (*,<*,< l»),  k(u)(»)  -  x»  f*#  If)***  ' 


where  k  is  defined  in  (5.6).  We  identify  the  variables  appearing  in  (8.17) 


with  those  of  (7.8): 

1  3_u  ^ 

(8.29)  a  (u(z),  (z),k(u)(z),z) 

AS 

2  3u  ^ 

ft  ( u ( z ) ,  (£>*]S<H:><5>»*> 

_  3u  ^ 

a  <u(z),  ~  (z),k(u)(z),z) 


A*  A#  A# 

+h(A)(z),z)  , 

A#  A#  AW 

aJ  /V  A/ 


31- 


^(ulz),  +  h<A>(£>'5)  » 

3u  M 

82(u(z) ,  (j5>#k(u)  (jg),z)  “  o(s)  , 

8(u{z),  (£)  »k(u)  (£)  >i)  ”  0  ' 

Y^u,*)  “  -Y(X(z),e(z),*(z),z)  , 

Y2(u»z)  *  -5(X(z),e(z),$(z),z)  , 

Y3(U/Z)  *  -5(X(z),*(z),*(z),z)  ♦ 

We  identify  w_  with  (X*,<p#,4>#).  Note  that  hypothesis  (8.22)  is  ensured  by  the 
restricted  strong  ellipticity  condition.  We  then  have 

A  A  A  A  M 

8.30.  Theorem.  Let  g»j§,fe»P/Y  ,«#8  satisfy  the  hypotheses  of  Theorem  8.20  with 
the  identifications  (8.28)  and  (8.29).  Then  (7.8)  with  x#  *  5  ia  satisfied  for 
all  (X#,<p#,x#)  in  V. 

The  question  of  regularity  of  solutions  for  the  types  of  systems  described 
in  this  section  remains  open.  Giaquinta  (1983)  gives  partial  regularity  results 
for  more  restricted  systems.  However,  it  is  by  no  means  clear  how  much 
regularity  is  physically  reasonable  for  either  of  the  more  general  types  of 
problems  presented  here.  Ball  (1982)  suggests  that  discontinuous  solutions  of 
problems  such  as  those  treated  in  (8.7)  can  be  used  to  model  rupture  of  solid 
bodies.  In  addition,  we  suggest  above  that  operators  such  as  (8.13)  can  be  used 
in  constitutive  equations  to  model  the  nonlocal  behavior  of  ferromagnetic 
materials,  and  the  physical  evidence  of  so-called  "domain  structures"  (cf.  Tebble 
(1969))  suggests  that  highly  discontinuous  magnetic  fields  are  to  be  expected 
from  a  good  model  of  such  materials. 


for  doing  so  ar«  virtually  identical  to  those  used  by  Antman  (1983).  The 
deformations  defined  by  (9.4a-c)  constitute  "family  2"  of  Truesdell  6  Moll  (1965 
Sec.  59).  The  other  functions  of  (9.4)  are  so  specified  as  to  ensure  that  our 
final  problem  consists  of  ordinary  functional  differential  equations. 

The  chain  rule  implies  that 

~  3X  3X  3x 

(9.6a)  £U(x))  -  (*(x))  “  3^  (x)  •  ^  (s(x) ) 

a*  —  m 

“  [w^ (s)e^ (jc)  +  w1  (s)w^(a)g2(x)  +  < s ) «3 3 is -| < 0 ) 

+  8’1t022W1(s)22(i>  +  °32S33~2(0>  +  [°23W1  (b)22(*)  +  a33S3l2S3  ' 
(9.6b)  3^£^x^  “  fg  *£****  *  w^( s)k^ ( 9)  , 

(9.6c)  e(*(x))  “  f|  (2<x))  -  v^(s)jc1(9)  +  +  053X3  * 

( 9 . 6d )  (z(x))  -  w^(s)k1(8)  +  a62s“1k2(0)  +  . 

The  representation  (9.6a)  reduces  (3.2)  to  the  requirement  that 

(9.7a)  (a22a33  *"  ^23^32)  >  0  a.e. 

Since  Wj  is  a  radial  distance,  me  require  that 
(9.7b)  w^s)  >0  for  s  e  [a,1]  , 

whence  (9.7a)  reduces  to 

(9.7c)  ^a22a33  “  a23a32*w1  *  0  a,e* 

Por  simplicity  we  require  that 

(9.7d)  w^  >  0  a.e.j 

the  opposite  case  corresponds  to  an  eversion  (cf.  Antman  (1979))  and  provides  no 
further  technical  difficulties. 

Mote  that  the  components  of  (9.6),  (9.4d,f)  with  respect  to  the  indicated 
base  vectors  and  dyads  are  independent  of  6  and  z.  (It  is  easy  to  show  that 
(9.4a)  and  (9.4e)  are  the  most  general  forms  whose  gradients  have  this 
property.)  We  denote  the  ordered  set  of  the  components  of  _T_  corresponding  to 
(9.4),  (9.6)  by  the  single  symbol 


(9.8) 


y{ b)  m  (Y^(s),.«.,Y^g(s)) 

-  (w^ ( s ) » w^ ( s )w3 (s ) , w3 ( s ) , 032®  ^w^(®)»fli22®  \®23*1^®)'®33  * 
w^(s) *w4(a) ,w£(s) »a52®  1 'a53'h1^ )  . 

We  define  the  physical  components  of  the  dependent  constitutive  variables  by 
(9.9a,b)  Ti;J<I/®>  5  Si<2i>*S<I.**>-JSj<e> '  S^X'®)  5  M'j*0)  ' 

A  A  A  A 

(9.9c,d)  d^XtS)  5  bj<I»®>  5  b,£»*>*kj{X'®)  ' 

(9.9e)  j^X'*)  5 

when  V_  has  the  form  corresponding  to  (9.4)  and  (9.6),  assuming  that  the 

A 

constitutive  functions  T,  etc.,  are  such  that  these  constitutive  functions  for 
the  physical  components  depend  only  on  X  and  s.  These  representations  are 

A 

valid  when  the  constitutive  functions  T,  etc.,  are  hemitropic  and  depend  on  z 
only  through  s.  They  are  also  valid  for  special  forms  of  aeolotropy.  It  then 

A 

follows  that  3$/3z  depends  only  on  X  end  s.  We  also  assume  that  a 
depends  only  on  a. 

We  now  obtain  an  alternative  representation  for  h  in  terms  of  3i)>/3z 
directly  for  the  semi-inverse  problem;  the  specialisation  of  the  results  of 
Section  5  does  not  yield  the  new  representation.  Substituting  (9.4f)  and  (9.9e) 
into  (3.10)  we  get 

(9.10)  0  -  j^Xt®)*®)  ' 

(9.11)  Ishjts)]'  -  sj3(x(s)»«)f  h|(s)  -  -jjtX*®)'®)  • 

Thus  h  must  have  the  form 

A# 

1  r1  * 

(9. 12a, b)  h^s)  -  w£(s),  h2(s)  -  s  1  (og2  -  /  tj3(X(t) ,t)dt]  , 

s 

r1  “ 

h3(s)  -  og3  +  /  j2(x<t),t)dt  . 


(9.12c) 


Therefore  h  can  be  written  as  the  sura  of  the  gradient  of  (9.4g)  and  an  integral 
operator  (cf.  (5.2)). 

Condition  (9.10)  may  be  regarded  either  as  a  restriction  on  x  or  else  as 

A 

being  identically  satisfied  by  virtue  of  choosing  the  constitutive  function 
to  be  the  zero  function,  in  which  case  the  material  is  incapable  of  conducting 
electricity  in  the  radial  direction.  In  the  former  case  we  assume  that  (9.10) 
can  be  uniquely  solved  for  w£  in  terras  of  the  other  elements  of  x* 

(9.13)  w£(s)  -  e1(x”(8)#s)  , 

where  x”  stands  for  all  the  components  of  x  except  w£.  A  sufficient 

A 

condition  for  (9.13)  to  be  equivalent  to  (9.10)  is  that  w£  »  be 

strictly  increasing  and  assume  both  negative  and  positive  values.  That  this 
function  be  strictly  increasing  is  ensured  by  the  strict  raonotonicity  of 

A 

e  **■  %(?_'£)•  Equations  (9.10)  and  (9.13)  are  also  equivalent  in  the  important 

A 

special  case  that  j-|(X»s)  has  the  same  sign  as  w^,  which  occurs,  e.g.,  if 

A 

;J(£,z)  *  J(JL*i£)£  where  J  is  a  positive-valued  scalar  function.  In  this  case 
(9.13)  reduces  to  Wg  *  0. 

Let 

( 9 • 14)  w  *  (w^,..« ,Wg ) ,  X  *  (o22,..., *63)  • 

A  A 

By  controlling  the  dependence  of  j2  on  j3  and  h2  and  h3  we  can  imitate 
the  development  of  Section  5  to  show  that  (9.12b,c)  can  be  uniquely  solved  for 
h2  and  h3  in  terms  of  the  other  variables.  Thus  we  can  replace  these 
equations  with 

(9.15a)  h2<s)  *  s-1ta62  +  ,s)^  ' 

(9.15b)  h3(s)  =  a63  +  • 

Alternatively  we  may  observe  that  (9.12b,c)  is  equivalent  to  an  initial 
value  problem  for  h2,h3.  If  we  assume  that  there  is  a  number  p  >  1  such 

A  A 

that  j2(£(*)»*)  and  j  3  [yji  * ) ,  •  )  are  integrable  on  [a,1]  when  xe  Lp(ta,1]) 


and  that  thara  is  a  number  such  that 

I  |3ja(X(a),s)/3h6|  <  K(w,a.) 
of8-1 

whan  X  €  lp(ta,1]),  than  tha  standard  theory  of  ordinary  differential  equations 
(cf.  Bala  (1969,  Sacs.  1.5,  1.6))  implies  that  (9.12b,c)  has  a  unique  absolutely 
continuous  solution  on  [a,1],  which  we  can  represent  by  (9.15a,b).  Bote  that 
this  result  does  not  require  restrictions  like  those  of  Theorem  5.3  on  the  sise 

A  A 

of  3;J/3h  and  on  the  size  of  the  domain.  The  Arzela-Ascoli  Theorem  implies  that 
W^((a,1))  a  w_(*)  »*■  <2 ( w( * ) ^ •)»«,*•)  «  C°([a,1J)  are  compact  (when  this 
construction  of  h2  and  h3  is  used). 

We  henceforth  assume  that  the  representation  (9.15)  is  valid  and  that  k2 
and  have  this  compactness  property. 

We  are  now  ready  to  mite  down  the  governing  equations  for  our  semi-inverse 
problem  when  the  only  body  force,  the  tjorentz  force,  is  absorbed  into  the 
effective  stress  and  when  the  only  heat  source  is  that  of  (3.6),  due  to  Joule 
heating.  Let  l,e.  be  the  axial  vector  corresponding  to  L.  Let  us  set 
(9.16)  _£  *  ( £ ^ , • • • , ) »  X  m  Ihj »n2»0 ,n^ ,hg,0 ) 

with 

A  A  A  A  A  A  A 

(9.17a)  C(w*  ,w,a,v(»  ),s)  =  (T^,w^T2^,T3j,q^,d^,b^)  , 

(9.17b)  n^(w*  ,w ,u ,v( *  ) , s )  -  T21w2  +  <*22*  T22  ^  **23^23  1 

(9.17c)  T)2(w' 1 f2.‘ 'JZJ I*  -  *3  • 

(9. 17d)  n4(w' ,w>,X#X(*)fS)  =  a52B~^2  *  *53^3  • 

(9. 17e)  h5(s)  -  -a  , 

where  the  arguments  of  the  constitutive  functions  appearing  on  the  right  sides  of 
(9.17a-d)  are  X»*  *nd  with  ®very  h1 ,h2,h3  appearing  on  these  right  sides 
replaced  by  w£,  a^s”1  +  1 s),  ag3  +  *3(v( •),«_, s)  respectively.  Wote 

A 

the  definition  of  $2*  In  line  with  the  remark  following  (3.12)  there  is  no  loss 


of  physical  content  in  taking  A 2  “  *3  “  0.  We  do  so  because  it  simplifies  the 
ensuing  analysis.  Then  by  using  the  componential  form  of  (3.12),  we  reduce  the 
governing  equation  (5.13)-(5. 16)  to  the  following  system  of  ordinary-functional 
differential  equations  for  w_,cu 

(9.18)  (s|r  -  sn  -  2 

A  A 

where  the  arguments  of  2  an(*  U.  are  w.'  • ) , s .  (We  have  introduced  our 

constitutive  functions  in  (9.17)  with  the  argument  so  as  to  avoid 

confusion  in  Section  10  when  we  take  certain  partial  derivatives  of  these 
functions . ) 

If  (9.10)  is  equivalent  to  (9.13),  then  w£  is  completely  determined  by  the 
other  components  of  w_  and  ct,  which  can  be  found  from  the  remaining  equations 
and  side  conditions.  We  accordingly  discard  the  fifth  equation  of  (9.18),  which 
is 

(9.19)  (sdj)'  -  so  . 

We  regard  this  equation  as  determining  the  o  necessary  to  maintain  the  semi¬ 
inverse  state  (9.4).  This  interpretation  of  (9.19)  smells  fishy,  but  is  in  fact 
quite  reasonable:  Consider,  e.g. ,  constitutive  equations  of  the  form 

“  Je^,  d^  —  De^ 

where  J  and  D  are  positive-valued  scalar  functions.  Then  (9.13)  and  (9.19) 
require  that  a  ■  0.  Thus  when  (9.10)  is  equivalent  to  (9.13),  we  shall  simply 
ignore  (9.19),  regarding  (9.18)  as  the  suitably  truncated  system.  We  shall 
comment  on  boundary  conditions  below. 

A 

If  j  1  is  the  zero  function,  then  we  need  take  no  action  with  respect  to 

(9.19) . 

We  now  specify  boundary  conditions.  Our  prescription  is  compatible  with  the 
formalism  of  Section  7.  On  the  cylindrical  face  s  ■*  1  of  3B  we  either  fix 


the  outer  radius: 


(9.20a) 


w,(1)  -  5,(1) 


where  w,(1)  is  a  given  positive  number  or  else  we  prescribe  the  traction: 
(9.20b)  ? !<_£*<  1),.w(1>, a., w.f),1>  -5,(1) 

where  5,(1)  is  a  given  number.  (More  generally,  we  could  replace  5,0)  with 
5 , ( w( 1 ) ,a)  where  the  new  5,  is  a  prescribed  function.  Since  only  minor 
technical  difficulties  are  introduced  by  such  a  replacement  in  this  and  other 
such  Neumann  conditions,  we  do  not  bother  to  pursue  such  generality.)  We  fix  the 
deformation  to  within  a  rigid  displacement  by  setting 


(9.21) 


w,(1)  -  0  , 


(9.22) 


w3(1)  -  0 


On  this  face  we  either  prescribe  the  temperature: 


(9.23a) 


w4d)  -  w4(D 


where  w4(1)  is  a  given  number  or  we  prescribe  the  heat  flux: 


(9.23b) 


54(w'(1),i(1),o,i(*),1)  -  54(D 


where  54(1)  is  a  given  number.  Finally  we  fix  the  data  of  the  potentials  <p 


and  <|»  by  taking 


(9.24) 


(9.25) 


Wed)  -  0  , 


w6(1)  -  0  . 


On  the  cylindrical  face  s  -  a  we  prescribe  alternative  boundary  conditions 
expressed  in  an  analogous  notation: 

A  mm 

(9.26 a,b)  w^fa)  -  w^{a)  or  5^(j»' (a),j»(a),£,wj{  •  ),a)  -  S^fa)  for  i  -  1,...,6 

where  w^(a),  i  -  1,...,6  and  5j,(a),  i  -  1,3, 4, 5, 6  are  given  constants  and 


‘.V.v. 


v*  •/> 

**  K«r , 


•>?>> 

\V.V?V 
*  V  V 

m 

CVTs.'-V" 


where 


(9.26c) 


52(a)  =  w,(a)T 


with  t  a  given  constant.  The  form  of  52(a)  reflects  its  definition  and  the 
fact  that  it  is  a  torque.  In  conformity  with  the  condition  that  wj  >  0,  we 
require  that  w, (1)  >  w, (a)  when  both  these  numbers  are  prescribed.  Note  that 


\ 

r  *  * «« 


-  v  v  vv.-.'  >  v/v 


i  V  •  *  •  ►  *  j 


•  »  1  «  s  .  *  _  > 

* • *  «  *  »  *  ’ 


(9.13)  implies  that 


(9.27)  w5(1)  -  w5(a)  -/  e1(x”(®)#s)ds 

a 

so  we  are  not  free  to  prescribe  both  w5(1)  and  wg(a)  when  (9.13)  holds. 

To  avoid  dealing  with  the  minor  technical  difficulties  that  can  arise  when 
all  boundary  conditions  are  of  the  Neumann  type,  we  assume  that  the  temperature 
w4  is  prescribed  on  at  least  one  of  the  faces  s  ■  a  and  s  *  1.  He  need  make 
no  such  provision  for  the  variable  Wj  because  the  growth  conditions  we  shall 
impose  on  our  constitutive  functions  preclude  any  trouble  with  coercivity 
ultimately  due  to  Neumann  conditions.  If  8  is  an  entire  tube,  then  9  «  *.  if 
we  require  that  y,  and  $  be  continuous,  then 

(9.28)  C»22  *  <*32  "  9#  ^S2  **  9  • 

We  obtain  various  kinds  of  dislocations  by  suspending  (9.28a,b).  If  8  is  a 
sector  of  a  tube,  i.e.,  if  0  <  v,  or  if  8  is  a  slit  tube  i.e.,  if  0  -  v  but 
with  the  faces  0  *  -it  and  0  *  ir  not  identified,  then  we  can  prescribe  certain 
degenerate  boundary  conditions  on  the  faces  0  *  ±0.  We  likewise  prescribe  such 
conditions  on  z  *  ±Z. 

We  adopt  the  following  alternative  conditions  for  the  faces  0  *  ±0: 
(9.29a,b)  a22  =  a ^  or 

A22^— 1 ■ J  5  /  w^(s)T22(]L(s) ,s)ds  -  ^22^1^  5  /  Wi(s)*22*s)ds 

“a  a 

where  a22  is  a  given  number  and  A22  a  9i-ven  functional  of  having  the 

A 

indicated  form.  In  the  argument  J{b)  of  T22*  h.j(s)  is  replaced  by  (9.12a) 

A 

and  h2(s),  h3(s)  by  (9.15).  It  is  easy  to  see  that  "ZA22^— ^  *-s  the 

resultant  effective  torque  about  e.  on  the  material  face  0  -  -0  needed  to 


maintain  the  stata  (9.4)  (cf.  Antraan  (1983)).  Similarly  we  prescribe 


(9.30a, b)  <*32  "  *^32  (v .a ]  =  f  s)de  ■  # 

a 

1 

(9.31a, b)  ®52  m  ^52  ^52 ~  J  <lj(Y( s)  ,s)da  ■  » 

a 

1 

(9. 32a, b)  a^2  "  ^62  A^2  [w.8 ]  a  J  bj ( Y  ( s )  ,s) ds  *  Xg2  • 

a 

A 

Here  X  ha*  the  form  just  described.  -ZA32  [w^aj  i8  the  effective  resultant 

force  in  the  e^- direction  on  the  material  face  0  •  -0.  ^32 '^52 »^62  are  3UBt 


numbers. 


For  the  faces  z  *  ±Z,  we  likewise  prescribe 


(9.33a,b)  o23  -  o23  or  A^lW/O.]  s  /  sw1(s)T23(x(s),s)da  -  *23tw15  ' 


(9.34a, b)  <*33  ■  °J3  or  A^ [w,a]  =  /  sT33(X(b) ,s)ds  *  S33  , 


(9.35a,b)  u53  -  S53  or  A^lW/O]  =  /  sd3(x( ■  )  »s)ds  -  Xjj  , 


( 9. 36a, b)  ag3  -  u63  or  Ag3t*/CL]  =  /  sb3(X(B)»*)<*s  - 


-QAjsfWjOj  is  the  resultant  effective  torque  about  e3  and  -0A33  is  the 

resultant  effective  force  in  the  e, -direction  on  the  face  z  ■  -Z. 

~  3 


10.  Consequences  of  the  Strong  Ellipticity  Condition 


A  A 

In  this  section  £  and  have  the  arguments  listed  in  (9.17).  Thus  the 

A 

derivative  of  h 1  with  respect  to  w1 ,  say,  is  a  pure  partial  derivative!  no 
differentiation  with  respect  to  v(«)  is  required. 

A  A 

Since  3T/3wJ  -  {3£/3p) se^k^ ,  etc.,  definition  (6.1)  and  the  strong 
ellipticity  condition  imply  that 

(10.1a)  v.(3i/3w')*v.  -  «>((v1e1  +  +  ^^I'^-l^rVl*  >  0 

for  all  v.  s  (v . , .  • .  ,Vg )  ^  _0_  when  w1  >  0.  Slightly  abusing  the  notation,  we 

likewise  obtain  A 

3£ 

(10.1b)  v«  t—  »v  >  0  W  +  0  when  w«  >  0  where 

X.  =  (w^wjw^w^w^w^w^)  . 

Next  we  observe  that 

2  A  A 

( '10.2)  wjcl  ”  a* h1  ”  a»T*c 

with  3  “  §2,  c  “  wySl  +  (a22/s)Jj2  +  023^3#  80  that  the  8tron9  ellipticity 

condition  implies  that 

A  A  A 

(10.3)  3n ./3w.  =  acs (3t/3w  )  =  (ac) : ( 3T/3F) sac  >  0  . 

I  I  ASA#  A/  |  ASA^  A/  A#  ASA# 

(Note  that  3  and  c  are  not  independent  of  a*F* c/| c | ) . ) 

Suppose  that  a32  end  a33  are  prescribed.  Set 

(10.4)  U  ■  *33° 22  "  ®32a23'  V  *  ®32a22  +  ®33a23  ' 

(10.5)  M  *  ®33T22  ”  ®32sT23'  N  *  ®32T22  +  ®33sT23  * 

We  solve  (10.4)  for  a22  and  a23  in  terms  of  y  and  v  and  substitute  the 

A  A 

resulting  expressions  into  the  arguments  of  H  and  N.  The  strong  ellipticity 
condition  then  implies  that 

ASM/Sy  3M/3v  \ 

(10.6)  I  A  ,  ]  is  positive-definite  . 

'‘SM/Sy  3n/3v/ 

The  combination  of  (10.1b)  and  Hypothesis  6.6  supports  a  global  implicit 
function  theorem  (based  on  degree  theory)  that  ensures  that  the  function 
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11.  Growth  Conditions  and  Function  Spaces 


We  introduce  some  notation  to  be  used  in  the  rest  of  this  paper.  I«et 


(11.1) 


l22  -  s  1w1^22'  ^23  5  w1*23'  ^32  “  *  1*32'  ^33  "  *33' 

A  A  A  A  —  jA  A  A 

C52  5  8  d2»  553  ”  d3'  C62  5  8  b2'  *63  "  b3  ' 

i*£L#  =  E  ^yva#yv  ' 


the  summation  being  taken  over  y  =*  2,3,5,61  V  -  2,3.  I<et 


(11.2) 


u  s  (w,a)  . 


We  set 


(11.3) 


a  A  A  A 

<m(w_)  ,Wjl ) >  =  J  (C«wi  +  n»w>  +  C_*a, ) ads 


a£(a).w#(a)  -£(1)^(1)  -  l  A  [w,]#  , 


where  the  arguments  of  5_,n_,C.  are  vf1  ,<»,<*.,  w( •  ) , s .  Observe  that 


(11.4) 


/  £»a^sds  -  l  Ryvl£i.I<x#yv  * 
a  yv 


The  weak  form  of  the  boundary  value  problem  of  Section  9  is 


(11.5) 


<m(u)  ,<jJj>  =*  0  "for  all"  . 


In  the  next  section,  we  give  precise  interpretations  to  relatives  of  (11.5). 

We  pose  the  basic  growth  conditions  in  terras  of  a  scalar  function  W,  which 
might  be  interpreted  as  a  sort  of  stored  energy  function.  It  allows  us  to 
replace  standard  Lp-spaces  by  related  spaces  better  equipped  to  handle  possible 
aeolotropy. 

11.6.  Hypothesis.  There  are  numbers  q  >  p  >  1,  c^  >  0,  C2  >  0,  c^  >  0,  c^  >  0 
and  a  function  R15  3  Wj.,*)  e  *  having  the  following  properties: 


^■-*.*.*  - .  •  .  ■  -  '  -  -  .  -  .•  .  ’  <_■  ...  .  •  .  '  ■■  -  •  -  *  .  *  c  -  *  V*  O 


mV 


(11.7) 


(11.8) 


W(«,s>  la  strictly  convex  , 


W(  0. s)  -  0  , 


(11.9)  w(  * ,s)  la  invariant  under  the  chan 


(11.10) 


cllx!P  -  °2  *  "<X»«>  <  c3lxlq  +  e4 


Let  w  be  in  affine  function  satisfying  whatever  Dirichlet  conditions  from 


(9.20)-(9.26)  are  prescribed.  Let  a  be  a  vector  of  the  form  a  of  (9.14)  with 


its  entries  taken  to  be  a......  whenever  these  numbers  are  prescribed  in 


Let  Y_  be  generated  from  and  b£  (9.8),  (9.12),  (9.15).  Then 

(11.11)  |mw'  -  w’)  +  £•  (w  -  w>  +!*(»  -  & )  >  w(x»s)  -  w(£,s) 

when  the  arguments  of  the  constitutive  functions  on  the  left  side  of  (11.11)  are 


w * , * w( * ) » s  and  with  X  and  i.  expressed  in  terms  of  these  variables  by 
(9.8),  (9.12a),  (9.15). 

Remarks.  Condition  (11.9)  is  a  sort  of  isotropy  condition.  Its  provenance  is 
described  by  Antman  (1983,  eg.  (7.7)).  The  mechanical  terms  from  (11.11) 
correspond  to  a  certain  stress  power.  This  issue  is  likewise  treated  at  great 
length  by  Antman  (1983). 

We  now  introduce  function  spaces  naturally  associated  with  W.  Let  n  be  a 
positive  integer.  We  set 

1 

(11.12)  G  s  {x  :  /  sW(x(s>i*0ds  <  "}  » 

a 

(11.13)  E  =  { (W/°L>  «X  €  G)  , 

(11.14)  En  *  ( ( W/«.)  €  E  t  nw}  >  1  a.e.,  nw^s)  >  s,  n(a22a33  “  a23°32)  >  ^  • 

(11.15)  A  5  {(W/CO  e  E,  w^  >  0,  w,  >  0  a.e.,  a22°33  “  °23032  >  °' 

a  fixed  subset  of  {w(a) ,_w( 1 ) ,(*}  is  prescribed  as  in  Sec.  9}  , 

(11.16)  An  5  A  n  En  * 

A  is  the  set  of  admissible  functions.  Conditions  (11.7)— (11.10)  ensure  that 


W(*,s)  satisfies  the  A, -condition  of  Orlicz  space  theory,  whence  it  follows 


that  6  is  a  reflexive  separable  Banach  space  satisfying 

(11.17)  Lq((a,1)>  C  G  C  I^<(a,1)) 

(cf.  Krasnosel'skii  &  Rutitskii  (1958)).  Since  some  components  of  X  are 
products  of  components  of  (w,a),  neither  E  nor  En  is  a  Banach  space.  It 
easy,  however,  to  construct  a  suitable  Banach  space  for  (w_,a):  Let  ^w,a_)  e  E 
let  be  defined  by  (9.8)  with  (w2, .  •  •  *Wg#&.)  replacing  (w2, .  •  *  ,Wg,a_) ,  an 

let  X  be  defined  by  (9.8)  with  w^  replacing  w^.  Define 

1  1 

(11.18)  Vn  =  { (^J/2D *  /  sW(X# (■)  ,s)ds  <  ®,  /  «W(X  (s),s)ds  <  •}  . 

a  a 


Proposition.  /_  is  a  reflexive  separable  Banach  space. 


and  A 


are  closed  subsets  of  y„. 


n  is  sufficiently  large.) 


The  proof  of  this  result  is  identical  to  that  of  Antman  (1983,  Prop.  7.25.) 
He  now  refine  (11.11): 

11.20.  Hypothesis.  There  are  positive  constants  o5,c6,c7,c8,e ,  depending  on 
wrw2,  such  that 

(11.21)  (wj  -  w1’)i1  +  (wt  -  w1)n1  >  c5(|w*jP  +  Iw^sP)  -  c6(1  +  lxlP_e>  • 

A  A 

where  the  arguments  of  are  those  listed  in  (9.17). 

The  preceding  hypotheses  ensure  that  the  material  is  not  too  weak;  the 
following  hypothesis  ensures  that  it  is  not  too  strong. 

thesis.  Let  the  constitutive  function  introduced  in  (9.9)  depend 


only  on  X*8,  Let  a(x)  and  a(jc)  be  vectors  that  are  linear  combinations  of 
£y&2'&3  with  coefficients  depending  only  on  s  and  let  g(]£)  and  c(x)  be 


with  coefficients  dependin 


only  on  s.  Let  a*  F*  c  *•  a*T*c,  g*a  w  g*a,  e*a  »  d*a,  h*a  »•  b»a  be  strictly 
increasing  and  let  t*(ac)  be  constant  when  g  has  the  form  (9.6).  Then  there 
are  continuous  functions  (£.,s)  *  T+(£_,^  ,  g+(I\z),  £+(£#,g)»  b+(T,z),  P+(T,z) 


M  ^  M  ^ 

with  a(x)»T  (£.(*)  »*)*c(x) ,  j3  (Hje)  ,*)»a(x>/  ate. ,  depending  only  on  s  when 
T  has  the  form  corresponding  to  (9.4)  and  (9.6)  euch  that 


(11.23) 


-v  +  ~  if  3p<*s>  +  *  * 

Si  * 

a»T(r,E)»c  if  9G(ac)  -  t  , 

*«v  —  /v  m  '  /saw 

+  I  *  ^  i  A  |  +  <w  |  A  a#  | 

a  #a  >  a*a  #  d  »a  >  d*a  ,  b  *a  >  b*a  , 

r#  a#  *w  '  m  <v  <«w  iv1  <v  «w  >  /w  «w  1 


p+  >  lk«i 


whan  £_  has  the  form  corresponding  to  (9.4)  and  (9.6).  Let  £.+  ,n>+,£_+  2J£. 


by  (9.17)  and  (11.1).  lf_  oi,ui^  e  E,  then  s  *  ^‘w^Js)  +  n,+*iS#^8^  +  C+*°tf(s) 
where  the  arguments  of  5.+/H.+  'JL+  are  w,' (a) #w.(b) #®.*Jl(* ) »*»  Is  Integrable  on 
a, 11.  Moreover,  if  to  is  confined  to  a  subset  of  E  corresponding  to  a 


bounded  subset  of  6#  then  the  corresponding  (£r,TT_T’,£_’’ )  and  m(w)  generate 
elements  confined  to  a  bounded  subset  of  G*.  In  particular 


(11.24)  i1  <  E),  n,  <  n,,  m  <  vr  , 

|c2l  <  \i3\  <  |«|  <  N+ 

where  s  E^(w,' (s)  ,w(s)  ,o_,w(* ) » s)  # . . .  are  in  L^darl])  and  are  confined  to  a 


bounded  subset  of  L«  when  their  ar 


subset  of  G. 

A 

Condition  (11.24)  restricts  the  response  of  n^  (as  well  as  other 
functions)  in  tension.  We  now  formulate  an  hypothesis  to  control  its  behavior  in 

A 

compression.  It  furnishes  a  quantitative  statement  of  how  n  j  is  influenced 
more  by  changes  in  w^  than  by  changes  in  w}.  Let  the  function  with  values 

t  -1  » 

ni (s  w1 #w4,£,o,v(» ) ,s)  be  the  composite  function  obtained  from  n 1  by  using 
(10.7b)  to  replace  its  first  set  of  arguments  w.'  with  those  of  in 
(10.7b).  Let  f*  be  f^  with  its  arguments  in  the  same  order  as  those  of  n*. 

thesis.  Let  there  be  a  number  C  >  0.  an  octuple  a,  a  scalar¬ 


valued  function  w.,  a  function  E  with  values  in  R  ,  and  a  function  v 


with  (v,o)  «  A  such  that 


(11.26)  |oJ  <  C,  °22a33  ”  ®23a32  *  ^  C»  |^J  <  C  , 

£•)(▼_'  'Z/2WL**  >»•)  <  c*  Uj(i* »▼/«/▼(•>#■)  I  <  C  for  j  -  2, ...  ,6 
Then  there  ia  a  number  m  >  0  (depending  on  C)  such  that 
(11.27a)  (0/m)  a  u  fl(u,w^  ,£,a,vlt  ),s)  is  decreasing  , 

1  4  * 

(11 .27b)  (0,m)  3  u  **•  Tt^(u,w4,£,c^,v(»),s)  is  increasing. 

Moreover 

1  t  _i  8  ♦ 

(11.28)  limsup  /  Hjte  +s  /  f1  (e  #w4(t)  ,C.(t)  ,v(  •  ),t)dt, 

e+0  a  a 

w4(s)  ,£(s)  ),s)ds  -  -» 

for  each  fixed  x  €  (a,1]. 

The  motivation  for  this  hypothesis  is  given  by  Antman  (1983). 

Our  final  growth  condition  is 

11.29.  Hypothesis.  There  is  a  positive  constant  Cq  such  that 


12.  existence  of  Classical  Solutions 


In  this  section  we  prove  that  when  the  Strong  Ellipticity  Condition  and 
growth  conditions  hold,  then  a  certain  set  of  the  boundary  value  problems  posed 
in  Section  9  have  regular  solutions.  We  restrict  the  data  prescribed  in  the 
alternatives  (9.29) ,  (9.30),  (9.33),  (9.34)  to  be  one  of  the  following  nine  sets 

(12.1)  (® 22 ,a 23 ,a 32 33 ^  ' 

AAA 

<A22'a23'°32'a33)'  (°22 r A23 '°32  '°33  h '  (<*22 '°23 'A32 '°33  >  ' 

AAA  A  A 

(O22'0l23'a32#A33>'  (A22,A23'a32'°33) '  (A23'°23'A32'a33)  ' 

A  A  A  A 

*a22'A23'a32'A33* '  *a22'a23'  *42^33^ 

because  the  unprescribed  variables  from  (a22'a23'a32'a33*  are  then  confined  by 

(9.7)  to  an  open  half-line  or  open  half-plane.  It  then  follows  that  the 

corresponding  An  is  a  closed  convex  subset  of  l/n. 

Since  our  present  work  generalises  that  of  Antman  (1983),  we  emphasise  only 
those  aspects  that  are  novel.  His  work  may  be  consulted  for  motivations  and 
further  discussion  of  such  matters  as  the  data  omitted  in  (12.1).  Our 
presentation  also  tacitly  corrects  some  flaws  in  his  arguments. 

Our  basic  result  is 

12.2.  Theorem,  bet  the  monotonicity  conditions  (10.1),  (10.3),  (10.6)  hold. 

Let  the  Growth  Hypotheses  6.6,  6.10,  11.6,  11.20,  11.22,  11.25,  11.29  hold.  Let 
one  of  the  sets  of  data  of  (12.1)  be  prescribed.  Then  the  corresponding  boundary 
value  problems  of  Section  9  have  classical  solutions. 


Proof. 

Step  I.  Existence  of  a  solution  to  a  truncated  variational  inequality.  We  can 
write 


(12.3) 


<tn( u_)  ,u^>  =  <n_( ,to^ > 


(12.4) 


<JL<!£.1  ^2)/^>  =  /  {i((w1),(9),w2(B),ct2,w2(»),s)«wfl(8) 

a 

+  H(  (wj2)  '  (s)  ,1rf2(S)  ,<*_2,W2(*  >  >8^*Wf  (8>}b<*3 
+  a|_(a)*w#(a)  -  £(1)^(1) 

+  E^vt“2l  “  Suvtwi3>a#Wv 

where  the  last  term  is  summed  over  w  *  2, 3, 5, 6,  v  *  2,3. 

Since  An  is  a  closed  convex  subset  of  V n,  since  (10.1)  ensures  that  jn 
is  semi-monotone  on  An  in  the  sense  that 

(12.5)  Cn^1  ,<j1_2 )  -  n(m_2 ,m,2 )  ,u>^  -  oj_2>  >  0  Vw_ 1 , w_2  €  An  , 

and  since  ♦  k2(“.»*  )/k3(“.»*  )  are  compact  by  assumption,  we  can  use  Hypothesis 
11.22  to  show  that  m_  is  an  operator  of  the  "type  of  the  Calculus  of  Variations' 
(cf.  Lions  (1969))  from  An  to  the  dual  space  y*  of  Vn.  Thus  m^  is  pseudo¬ 
monotone  on  AR.  Hypothesis  11.6  ensures  that  m_  is  coercive.  A  theorem  of 
Brezis  (1968)  (cf.  Lions  (1969,  p.  297))  then  implies  that  for  n  sufficiently 
large  there  exists  an  €  An  satisfying  the  variational  inequality 

(12.6)  0  >  <m(o)_) ,(u„  -  in  >  Vti)  e  A„  . 

—  — n  —  —  n 

Step  II.  Bounds  on  Y^. 

Let  u  »  ■  (w^*®,)  where  (w ,a)  is  defined  in  Hypothesis  11.6.  Let 

and  correspond  to  and  oi^.  Then  (12.6)  and  (11.11)  imply  that 

1 

(12.7)  0  >  /  {£  • (wi  -  w’)  +  (w_  -  w)}sds 

a 

+  af (a)* (w  (a)  -  w(a)] 

Tl 

-  C(1)*  [w  (1)  -  w(a)] 

—  “T1  — 

+  ^  ~  \v  ^wn1  ^  ^any  v  “  ayv^ 

H  ,v 

1  1 

>  /  Wty^ts)  ,s) sds  -  /  W(y.(s)  ,s)sds  -  A 

a  a 

A 

where  C^ts)  =  5.(w^( s)  ,1^ s)  .a^w^* )  ,s) ,  etc.,  and 


J 


(12.8) 


1 

A=-/  ctaHw^U)  -  w5(s)]sda 

a 

+  £(1)‘ (w  (1)  -  w(D]  -  a£(a)*[w  (a)  -  w(a)] 

"  *Tj 

+  1  ^uv  ^Wn1^°nwv  “  auv*  * 

M/V 

In  the  following  development  we  let  C  represent  a  positive  constant 
independent  of  n,  which  can  always  be  estimated  in  terms  of  the  available 
data.  The  meaning  of  C  can  vary  with  each  appearance.  Note  that  (9.21) 
implies  that  wn2(1)  *  w2(1)  ■  0.  Conditions  (9.26c)  and  (9.21),  the  positivity 
of  w^1 ,  and  the  Holder  inequality  then  imply  that 

(12.9)  |a|2(a)[wn2(a)  -  w2(a)] | 

<  Cawn1(a) [|wn2(a)  -  wn2(1)|  +  1] 

1 

<  C  /  wn1(8)!w^2<s>!8d®  +  ^nl^ 

a 

<  Cl  w  ,w'  ,b  I  +  Cw  « ( 1 )  . 

nl  n2  p  **1 

In  this  way,  by  using  (9.21)-(9.26),  the  Holder  and  PoincarS  inequalities,  and 
the  estimates  (11.30)  we  find  that 

(12.10)  A  <  C{ 1  +  lY  ,L  I  +  Iw  ,w’ I }  . 

■ti  p  m  p 

(For  certain  sets  of  data  in  (12.1),  Iw  ,,W^I  <  ClY  ,L  I.)  Combining  (11.11) 

m  p  “Ti  p 

with  (12.7),  (12.9)  we  obtain 

(12.11)  lY  ,L  lP  <  C{1  +  Iy  ,L  I  +  Iw  ,,w'l}  . 

p  hi  p  nl  p 

m  a# 

Next  we  take  all  the  components  of  (£_  except  w^  to  equal  the  components 

»w  ^ 

of  &>n  and  we  take  w^  *  w^.  Then  (12.6)  and  (11.22)  likewise  yield 

(12.12)  Iw  4,wV  <  C{  1  +  lY  ,L  1  P”e  +  Iw  ^W1!}  . 

m  p  p  m  p 

Inequalities  (12.11)  and  (12.12)  imply  that 

(12. 13a, b)  lY  ,L  I  <  C,  Iw  ..W1!  <  C  , 

“Ti  p  m  p 
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whence  A  <  C.  It  follows  from  (12.7)  that 


(12.14) 


/  W(ln(s),s)sds  <  C 


Thus,  by  the  definition  of  the  norm  of  G  (by  duality  according  to  the  theory  of 
Orlicz  spaces),  we  obtain 

(12.15)  *1*'®'  *  c  * 

We  accordingly  get  corresponding  bounds  on  all  the  components  of  except 

wn2,  <*n22»  an23*  70  bound  these  variables  we  need  a  uniform  positive  lower  bound 


Step  III.  Integral  inequalities.  We  now  make  judicious  choices  for  in 
(12.6)  in  order  to  extract  useful  consequences  from  it. 


wn2' 

an22' 

for 

Wn1  * 

Step 

p  III. 

(12. 

6)  in 

If  l 

and 

set 

(12.16) 


wni(s)  for  a  <  s  <  x  -e  , 

w^ (s)  =  ‘  wn1(s)  +  [s  -  (x  -  e)]/e  for  x  -  e  <  s  <  x  , 

W_.,(s)  +  1  for  x  <  S  <  1  , 


(12.17)  w  *  (w<|(s),wn2,...,wn6),  o  -  . 

Then  w  *  (w,a)  €  when  w  and  a.  are  given  by  (12.16)  and  (12.17). 

Substituting  (12.16),  (12.17)  into  (12.6)  and  letting  e  +  0  we  obtain 


(12.18) 


x£nl(x)  >  C1(1)  -  /  snn1(s)ds 


for  almost  all  x  in  (a,1).  Since  the  right-hand  side  of  (12.18)  is  a 
continuous  function  of  x,  we  can  assume  that  (12.18)  holds  for  all  x  in 
(a,  1 ) . 


•  "  4  •  • 


.vy.  V V\  .'vN. 


**  fiT'  *v  •. 


\  \* 


“•****•  V  V  V  v  W  A 

-v  *v.\ *•' •> a\v£v v^\v>- 

.  *\  m*  w  .  _  V*  W  «  •» 


Zf  on  the  other  hand  w^(1)  is  prescribed  to  equal  w^(1)  by  (9.20a),  then 
we  require  a  more  delicate  construction.  Let  PR  be  the  set  of  all 
y  €  (a,1]  for  which 


(12.19) 


y 

lim  e  1  /  w'.(s)ds 
e+0  y-e 


exists  and  exceeds  1/n;  let  be  its  complement  in  [a,  1  ] .  (The  theory  of 

differentiation  ensures  that  (12.19)  exists  a.e.  On  p£,  wr1(s)  »  1/n  a.e. 

Below  we  show  that  the  Lebesgue  measure  of  approaches  0  as  n  ♦  «.  For 

now,  all  we  require  is  that  Pn  not  be  empty.) 

Let  us  choose  y  e  Pn,  x  t  (a,y)  and 

(12.20)  0  <  1  <  wRl(y)  "  wnl<x*  *  “  (y  ”  *)  • 

the  rightmost  term  of  which  is  positive  by  the  definition  of  Pn«  Since 

(12.21)  s  **•  wn1(s)  +  (y  -  s)/n  =  en(s) 

is  continuous,  since  Sn(y)  ■  wn.,(y)  >  wnj(y)  "  ®nd  since  »n(x)  <  wRl(y)  “  * 
by  (12.20),  the  intermediate  value  theorem  ensures  that  the  equation 

(12.22)  wni< t)  +  X  -  vn1(y)  -  (y  -  t)/n 

has  a  solution  t^X)  «  lx,y).  Since  sn  is  nowhere  decreasing  on  [x,y] ,  all 
solutions  of  (12.22)  lie  in  an  interval,  which  is  closed  because  *R  is 
continuous.  Since  t^X)  satisfies  (12.22)  and  since  y  e  Pn,  it  follows  that 
there  is  a  positive  number  6,  depending  on  wR^  and  y,  such  that 


(12.23)  ty  -  t^Xmn”1  +  9]  <  / 

P_ 


(X) 


W^tsjds  -  [y  -  tj^tXJl/n  +  X 


which  implies  that 


tn(X)  ♦  y  as  X  ♦  0  (for  fixed  n)  . 


(12.24) 


When  w1(1)  is  prescribed  to  equal  w.,{1),  we  let  y  e  PR,  x  «  (a,y), 
e  €  (0,x  -  a)  and  take 


(12.25) 


w.jts) 


wn1(s)  for  a  <  s  <  x  -  e  , 

wni(s)  +  [s  -  (x  -  e)]X/e  for  x-e<s<x, 
wn1(s)  +  X  for  x  <  s  <  tft(X)  » 


(.  wn1(y)  -  (y  -  s)/n  for  tR{X)  <  s  <  y  . 

We  define  «o  by  (12.25),  (12.17),  observing  that  u  e  AR.  We  substitute  this 

u_  into  (12.6),  let  e  +  0,  and  then  let  X  +  0  to  obtain 


(12.26) 


x£ni<x>  *  yCn1(y}  "  /  ®nn1(s)ds 


for  all  y  in  Pn  and  for  all  x  e  (a,y). 

By  the  simpler,  classical  version  of  the  process  leading  to  (12.18)  or 
(12.26)  we  likewise  obtain 


(12.27) 


xC_a(x)  -  Z-AD  -  /  snn..<s)ds,  j  ”  2,. ..,6  , 


for  almost  all  x  in  (a,1). 

To  be  specific  in  the  rest  of  our  analysis,  we  suppose  that  and  “33 

are  prescribed.  (Thus  we  can  exploit  (10.4)-(10.6).)  By  substituting 
»  *  (w  ,<*)  with 

(12.28a)  a  “  t®n22  +  ®33'  °n23  "  ®32'®32'®33* 

and  with 

(12.28b)  a.  -  (°n22  +  *®32'  ®n23  +  *®33'®32'®33* •  *  €  R 

into  (12.6)  we  obtain  from  the  arbitrariness  of  X  that 
(12.29a)  ®33A22 "  *32^3^^  >  ®33^22^wn1^  "  ®32^23^wn1^  • 

®32A22 +  ®33A23^^  “  ®32A22^wn1^  +  ®33A23^wn1^  * 


•/<*  -  •  •/  -  *  */«•  -  *  ./  «  A/  ,  t  «  •  .  ■  .  -  .  *  . 


•  •“  «rw  1 


PVV 


_ _ 

S§!» 

l my.* 

;r5 

*,j 

Step  IV.  Uniform  lower  bounds  for  Wn1'wn1'wn*  Inequality  (12.15)  enables  us  to 
use  the  second  part  of  Hypothesis  11.22  supported  by  (10.2),  (10.3)  to  deduce 
from  (12.18)  that 

(12.30)  *5ni(x)  > 

for  almost  all  x  in  (a,1],  and  from  (12.26b)  that 

(12.31)  x5n1(x)  >  tfnl^*  "  c 

for  almost  all  y  in  Pn  and  for  almost  all  x  in  (a,y).  Since  (12.30)  and 

(12.31)  yield  essential  lower  bounds,  we  regard  these  inequalities  as  holding  for 
all  such  x  and  y.  Since  h2  *  0  ■  n3,  and  since  (12.27)  consequently  implies 
that 

1 

(12.32)  (1  -  a)£nj(1)  *  /  s£nj(s)  ds  for  j  “  2,3  , 

ft 

we  obtain  from  (12.15),  (12.27),  and  Hypothesis  11.22  that 

(12.33)  |sCn2(s)l,  Is^ts)!  <  C  Vs  e  [a,1]  . 

Inequalities  (11.26),  (12.13b)  enable  us  to  deduce  from  (12.29)  that  there  is  a 
C  >  0  such  that 

(12.34a)  ®33a22^Ui^  "  ®32A23^£n^ 

=  /  wn1(s)M(J][n(s),s)ds  >  -C  , 

a 

(12.34b)  !®32*22^^  +  °33A23^n^l 

5  |/  wn1(s)N(j[n(s),s)ds|  <  C  . 

a 

The  analysis  in  the  rest  of  this  step  relies  critically  on  Hypothesis 
6.10.  He  identify  the  basis  {£T}  thusi 


(12.35) 


$1  *  JS)lSi  •  £ 2  m  ®3  *  $1^3*  ^4  "  ^2~1  ' 


~2(-C,33^2  ”  a3^S3) 


S2(“a32~2  +  0I33^3) 


Ss"  <.2;2.s2>1/2  '  <.2;2.;V'2 

33  32'  lB  32  33' 


S7  -  Sa&v  £b  ■  -8 3^2'  Sg  "  83JS3 


The  dual  basis  is  given  by 


(12.36) 


,2-  2  -  2. 1/2 
*  .  (S  *33  *  a32} 
m  -  2  -  2 
*«>32  *  “331 


e_(ak.  -  sale, ) 
~2  33~2  12«3 


32~3  ' 


,2-  2  -  2,1/2 

.*  (  32  +  33J  .t  .  .  -  .  . 

h - — - -- - S2«*„!S2  +  »«„£,>  ' 


”<S32  *  S33> 


Sx  m  £r  for  T  +  5,6  . 


(12.37) 


(12.38) 


JEsE-t  -  wi'  £^2  “  °'  ~‘®3  "  °'  ~*54  -  *1*2  ' 

F:F5  -  viW(s253^  +  5322>’V2'  £s56  -  *iV(825  32  +  5  \3)~'/2  , 

F:E  m  wii  FsE  m  8  U  ,  F:E_  ■  (X..  , 

~  ~7  3'  ~  ~8  32  ~  ~9  33 

T:E4  “  T . .  “  C4,  T:E  “  T4„,  TiE-  ”  T4_,  TiE.  «  T  , 

~  ~1  11  1  ~  ~2  12  ~  ~3  13  ~  ~4  21 


,  2-  2  -  2.1/2  ,  2-  2  .  -  2. 

(*  °33  *  a32>  *  *  •  (*  °32  *  a33>  * 

e  ^  a  Wj  T  t  ^  e  —  a  ”  > 

■,532  *  °33>  *‘a32  *  -33 » 


T:E_  *  T_„,  TiE„  *  T„,  T:E_  -  T__  . 

~  ~7  31  ~  ~8  32  ~  ~9  33 

(We  could  alternatively  take  E^  ■  e^c/ | c |  where  c  is  given  by  (10.2).) 

We  now  obtain  a  lower  bound  w*  for  wnl  that  is  positive  on  (a,1J. 
Suppose  for  the  sake  of  contradiction  that  wn1  have  no  such  lower  bound.  Then 
there  would  be  an  x  in  (a,1]  and  a  subsequence  such  that  *nj(x)  ♦  0. 

Then  wn1  *  0  uniformly  on  [a,x].  The  representation  of  wn1(x)  -  wn1(a)  as 
an  integral  of  w^  over  [a,x]  shows  that  Wn1  *  0  in  Lj(a,x),  whence 


m 


w 

v*" 

W  ’ 


ms 


has  a  further  subsequence  with  ♦  0  pointwise  a.e.  on  [a,x] .  zt  follows 
from  (12.15)  that  {1,5}  C  b,  {2, 3, 4, 6, 7, 8, 9}  C  e  for  almost  all  y  e  [a,x]. 
Condition  (12.34a)  ensures  that  no  alternative  of  (6.12b)  is  tenable  for  almost 
all  y  e  [a,x].  Hence  there  is  a  function  w*  such  that 


(12.39) 


w_.(s)  >  w«(s)  >0  Vs  e  (a,1]  and  Vn  . 


Note  that  we  can  define  w^  by 


(12.40) 


w«(s)  =  inf  w  . ( s )  . 


It  follows  that  w^  is  nowhere  decreasing,  for  if  x  <  y,  then 


(12.41) 


w?(y)  *  w!(x)  ■  inf  W  .(y)  -  inf  w n1(x) 


-  inf  wn1(y)  +  sup(-wn1(x)) 
>  inf(wn1(y)  -  wn1(x))  >  0 


since  wn1  is  nowhere  decreasing.  Xf  w^(a)  is  prescribed,  then 
w*(a)  ■  Wj(a)  >  0.  Otherwise,  we  have  yet  to  show  that  w*(a)  >  0. 

A  single  version  of  the  preceding  arguments  shows  that  there  is  a  C  >  0 


such  that 


(12.42) 


Un  >  VC 


We  now  confront  the  weakness  of  (12.26)  and  (12.31)  inhering  in  the 
membership  of  y  in  Pn,  which  conceivably  could  be  sparsely  distributed  over 
ta,1J. 

12.43.  Lemma.  Let  w„.(1)  ”  w. (1).  The  Le be ague  measure  of  Pf  approaches 


0  as_  n  +  •. 

Proof.  Were  the  conclusion  false,  there  would  be  a  C  >  0  and  a  subsequence 
{Mjj}  such  that  the  measure  of  P®  exceeds  c"1 ,  Condition  (12.15)  implies  that 
for  any  e  t  (0,  j  (1  «  a))  there  is  a  subset  of  [a  +  e,1]  of  measure 

1  -  a  -  2e  and  a  positive  number  Cg  such  that  the  absolute  value  of  each 
component  of  j^ts)  except  w^(s)  bounded  by  Ce  when  s  «  Q^,  Now  we 
choose  e  so  small  that  the  measure  of  P?  O  (^  =  R®  has  a  positive  lower  bound 


(independent  of  n).  We  fix  e.  Since  wn1  <  w(  1 ) ,  inequality  (11.24)  implies 


that 

(12.44) 


£  wn1Mnds  *  c  » 


so  that  (12.34a)  yields 

(12.45)  J  wn1Mnd8  >  * 

R 

n 

The  properties  of  Qg  and  the  alternatives  of  (6.12b)  imply  that  there  is  a 
sequence  «n  of  numbers  with  ♦  •  as  n  •  such  that  Mjj  <  -mn  on  R°. 

But  this  inequality  is  incompatible  with  (12.45)  if  the  measure  of  R°  has  a 
positive  lower  bound  independent  of  n.  □ 

We  next  obtain  a  lower  bound  K *  for  £n1  that  is  continuous  on  [a,1). 
This  bound  is  given  by  (12.30)  when  € 1 ( 1 )  is  prescribed.  We  prove  this  by 
showing  that  for  any  z  e  [a,1)  there  is  a  positive  real  number  h(z)  such  that 

(12.46)  ?n1(s)  >  “h(z>  Vs  e  [a,z]  . 

By  choosing  a  sequence  of  such  z's  approaching  1  we  obtain  a  sequence  of 
constant  lower  bounds  for  £n  whose  graphs  are  horizontal  line  segments  in  the 
( 8,Cni ) -plane.  By  joining  parts  of  these  segments  with  straight  lines  we  readily 
construct  a  lower  bound  continuous  on  the  half-open  interval  [a, 1 ) .  Suppose 
that  for  given  z  there  were  no  such  h(z).  Then  there  would  be  a  sequence 
xn  e  (a,zl  such  that  Sni(xn)  as  n  ♦  •.  Then  (12.31)  would  imply  that 

(12.47)  5n1(yn)  ♦  —  Vy„  €  PnO  (z,1]  , 

so  that  lea  for  (jSn(yn) •  Since  yn  is  bounded  by  (12.15),  (12.39), 
(12.42),  condition  (6.12a)  would  imply  that  M^y^  ♦  -*•»  which  is  impossible  by 
an  argument  like  that  centered  on  (12.45).  Thus  there  is  a  function  £* 
continuous  on  [a,D  such  that 

(12.48)  Cn^(s)  >  €*{s)  Vs  e  [a,1],  Vn  . 
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We  now  obtain  a  lower  bound  for  {w^} •  Iiat  _f_  be  defined  by  (10.7).  Then 

(12.48)  implies  that 

(12.49)  wn1(s)  “ 

>  f1(5*(s),Cn2(s),...,Cn6(s),wn(e),on#wn(»),s) 

>  min  f  ^ (C  ^  ( s) /C ^ ( s )  / . . . )  =  4*n(s)  . 
k<n 

\J»n  is  continuous  on  (a,1).  The  equation  (12.32),  the  bounds  (12.33),  and  their 

analogs  for  £n4,Sn5,£n6,  the  embedding  theorem,  the  bounds  on  un  and  vn, 

and  inequality  (12.39)  all  show  that  on  any  compact  subset  of  (a,1)  the 

sequence  (fn)  is  uniformly  bounded,  bounded  below  by  a  positive  function, 

* 

equicontinuous ,  and  decreasing.  The  Ar  zela-Ascoli  Theorem  implies  that  the  whole 
sequence  in  converges  uniformly  on  any  compact  subset  of  (a,1)  to  a 
continuous  limit  function  which  is  positive  on  (a,1).  (If  5^(1)  is 

prescribed,  then  is  positive  on  (a, 1]j  if  w^(a)  is  prescribed,  then 

♦  *  is  positive  on  [a,1).)  We  thus  have 

(12.50)  wn1(s)  >  • 

a 

Step  V.  Classical  solutions.  Let  0  <  e  <  j  (1  -  a).  Let  g  be  any  piecewise 
continuously  differentiable  function  with  g(s)  ■  0  for 
s  €  [a, a  +  e]  U  [1  -  e,1]  and  with  |g'|  <  y  <f>*.  Set 

(12.51)  *  9.  «  -  <V’,»1'-")  • 

Then  for  n  sufficiently  large,  Z  e  A„  (since  i|»*  has  a  positive  lower  bound 
on  [a  +  e,1  -  e]).  we  substitute  (12.51)  into  (12.6)  and  use  the  arbitrariness 
of  g  to  obtain  in  place  of  (12.18)  and  (12.26)  the  equality 

xcn1(x)  -  yCn1(y)  -  /  snn1(s)ds  Vx,y  €  [a  +  e,1  -  e]  . 

x 


(12.52) 


We  can  now  show  that  €nj(D  is  bounded  below*  Dots  that  ( 10.2) ,  (10.5) 
imply  that 

(12.54)  r»n1  "  wn1wA2«n2  +  B'^S  32  +  5  !3>"1^nMn  +  Vn>  * 

Hypothesis  11.22  or  estimate  (12.34b)  together  with  (12.39)  imply  that  { j  Nn  | } 

is  bounded  by  a  fixed  inte gratis  function  of  s  on  any  compact  subinterval  of 
(a,1J.  Hypothesis  11.22  and  estimate  (12.34b)  imply  the  same  for  |Mn| . 
Hypothesis  11.22  alone  implies  the  same  for  £n2«  Our  estimates  for  wn#vn  and 
our  representation  for  w^2  then  show  that  for  each  x  c  (a, 1]  there  is  a 
number  h(x)  such  that 

1 

(12.55)  |/  nn1(s)ds|  <  h(x)  . 

x 

Combining  this  estimate  with  the  limiting  equality 

(12.56)  x£  1  (x)  -  yC^y)  -  J  sn^sjds  for  a  <  x  <  y  <  1 

x 

corresponding  to  (12.26)  we  find  that  5^(1)  >  which  implies  that  £n1(1) 

is  bounded  below. 

To  show  that 

(12.57)  51(j»'(e),wj(s),£fw>(*  )#s)  ♦  ^(a)  as  s  ♦  a 

when  f 1 (a)  is  prescribed,  we  require  a  positive  lower  bound  for  wn1(a)  so 
that  we  aan  choose  an  2  with  w^ (a)  <  wnl(a)  for  large  enough  n.  Indeed, 
using  (12.56)  in  the  limiting  form  of  (12.6)  we  obtain 

(12.58)  0<  t$,(a)  -  51(a)][w1(a)  -  w^a)] 
so  that 

(12.59)  5,(a)  <  C1 (a)  . 
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He  can  now  use  (12.56)  to  conclude  that  (12.55)  holds  with  x  ■  a.  Thus 

(12.60)  Cn1(s)  <  *,<*>  +  C  <  C  . 

We  now  prove  the  existence  of  a  lower  bound  for  wr1 ( a) .  Since  f 1  is 
increasing  in  we  obtain  from  (12.59)  that 

(12.61)  *Al(s)  <  fJ(8’1wn1(8>'wn4(8)'c'5n2(8)'“-'5n6(8>'2n{,)'8)  * 

Now  suppose  for  the  sake  of  contradiction  that  there  is  a  subsequence  for  which 

wR(a)  ♦  0  as  n  ♦  ».  Estimate  (12.39)  ensures  that 

(12.62)  s^w^ts)  >  s_1w*(s)  >  w*(s)  >  a"1wn1(a) 

for  this  subsequence  for  n  sufficiently  large.  Property  (11.27a)  enables  us  to 
deduce  from  (12.60),  (12.61)  that 

(12.63)  wAl(s>  <  fI(a~1,,n1<a)*"n*<8)',’,> 

for  s  sufficiently  close  to  a  and  for  n  sufficiently  large.  Hence 

•  . 

(12.64)  wn1{8)  <  wn1(a)  +  /  »n1<a)  »%4(b) , . . .  )dt  . 

a 

Property  (11.27b)  then  yields 

s 

(12.65)  0n1(s)  <  n^(a_1wn1(a)  +  s  1  /  f^ (a_1wn1(a) ,wn4(t) , . . . )dt,wn4(s) , . . . )  . 

But  then  (12.55)  contradicts  (11.28).  Zt  follows  that  w^(a)  >  1/C  and  that  we 
can  consequently  choose  w^(a)  to  reverse  the  inequality  in  (12.58).  Hence 
(12.57)  holds.  The  demonstration  that  other  Neumann  conditions  at  a  are 
satisfied  is  routine.  These  results  complete  the  proof  of  Theorem  12.2. 
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